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SUMMARY 


Studies  of  tho  interface  properties  of  MIS  devices  made  with  photox"  si02 
an  BgCdTW  have  shorn  that  interface  trap  charge  (Qit)  is  an  important  contri¬ 
butor  to  net  interface  charge  In  soae  cases  it  can  dominate  the  fixed  charge 
resident  in  the  insulator.  This  indicates  that  the  net  interface  charge 
which  is  the  doe inant  surface  effect  on  passivated  detector  properties,  trill 

f  1 

be  dependent  on:  pre-Si02  deposition  surface  treatments;  water  content;  and 
surface  ferni  energy  level  ()  and  therefore  on  doping  density  and  temper¬ 
ature. 

The  implantation  task  has  demonstrated  an  ability  to  detect  atomic  den¬ 
sities  of  1H,  hi  and  *Hs  in  HgCdte  and  Cdie  above  101®  cm-3.  Cesium  beam  SIMS 
was  used,  and  range  parameters  have  been  determined.  This  will  be  an  impor¬ 
tant  diagnostic  tool  for  later  examination  of  channeled  Be  and  H  Isotopes. 

Angle  resolved  photoeaission  spectroscopy  has  been  used  to*  obtain  a 
detailed  description  of  the  bulk  HgCdte  electronic  band  structure.  The  elec¬ 
tronic  structure  of  the  alloy  at  once  accounts  for  the  desired  tunable  bandgap 
of  the  materials  as  well  as  many  materials  problems  in  growth  and  prepara¬ 
tion.  These  are  necessary  to  explain  the  surface  and  interface  properties 
affecting  passivation  as  well  as  the  LPE  heterojunction  interface.  Results 
are  also  presented  using  TEM  to  examine  this  interface.  Analysis  of  the 
imaged  defects  yields  data  an  stacking  fault  energies  and  dislocation  motion 
and  annealing.  Surface  properties  were  probed  directly  in  a  comparative  study 
of  activated  oxygen  uptake  on  the  clean  (110)  faces  of  CdTe,  HgTe,  and 
H9o.69Cd0.31T*'  *****  attempt  at  learning  more  about  the  mechanism  of 

surface  type  conversion  observed  upon  cleavage  by  performing  the  cleaving  at 
reduced  temperature.  Initial  data  regarding  the  kinetics  of  the  oxidation 
process  were  obtained  which  augment  information  on  the  chemistry  of  the  oxy- 
gen-BgCdTe  interaction. 

Bonding  studies  have  concentrated  on  the  calculation  of  important  alloy 
properties.  These  include  atonic  term  values,  electron  extraction  energies, 
and  many  body  excitation  energies  for  all  the  group  II,  III,  IV,  V  and  VI 
elements  which  are  important  for  tetrahedral  semiconductors.  Calculations  are 
also  in  progress  which  will  yield  the  bond  energies  and  bond  lengths  of  the 
"pure"  semiconductor  (binary)  compounds  and  their  modification  in  the  pseudo¬ 
binary  alloy. 
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•  interim  report  presents  results  obtained  in  the  second  six-month 


periodCJ  ending  December  31,  1983)  of  a  three-year  study  of  the  Interface 
region  of  BgCdVs.  The  effort  is  a  collaboration  of  groups  at  Santa  Barbara 
Research  Center,  Stanford  University  SRI  International,  and  Hughes  Research 
Laboratories.  The  technical  effort  is  divided  into  tasks  to  study:  1)  inter¬ 
face  properties  which  are  important  to  device  passivation;  2)  surface  effects 
and  bulk  defects  due  to  ion  implantation;  3)  properties  including  defects  of 
the  free  surface  and  defects  of  the  alloy  itself;  and  4)  theoretical  analysis 
of  bonding  in  random  substitutional  alloys  as  well  as  defect  formation. 


Results  of  these  tasks  are  presented  in  the  following  four  sections  of  this 


report. 
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Section  Z 
Interface  Study 


The  purpose  of  the  interface  task  is  to  determine  the  electronic  charac¬ 
teristics  of  passivated  interfaces  concentrating  on  properties  which  have  the 
greatest  effect  on  device  operation.  We  are  concerned  principally  with  opera¬ 
tion  of  photovoltaic  hetero-  or  homoj unctions  of  the  small  size  seen  in  large 
FPA  applications,  these  are  the  applications  where  surfaces  can  have  the 
dominant  effect.  As  discussed  in  the  previous  interim  report  the  principle 
impact  on  diode  RQA  is  due  to  interface  charge  which  can  produce  a  depletion 
layer  at  the  HgCdlte  surface.  Interface  charge  is  viewed  as  being  due  to  a 
combination  of  fixed  charge  (in  the  dielectric),  slow  interface  trap  charge, 
and  fast  interface  trap  charge,  depending  on  the  response  time  of  the  particu¬ 
lar  states  involved.  These  are  indicated  by  a  shift  in  flat  band  bias  a 
voltage  hysteresis  (seen  by  comparing  the  voltage  shift  of  the  accumulation  to 
inversion  sweep  with  the  return  sweep) ,  or  a  distortion  of  the  slope  of  the 
high-frequency  response  through  depletion. 

The  amounts  of  these  types  of  charge  vary  strongly  from  device  to  device 
for  what  are  nominally  the  same  interface.  For  the  case  of  Photox”  Si02  on 
HgCdTe,  the  Si02  is  altered  depending  on  (1)  pressures  and  mixtures  of  the 
reactant  gases,  deposition  rate,  substrate  temperatures;  (2)  the  details  of 
the  surface  preparation  (i.e.,  whether  a  polish  was  used  in  conjunction  with 
an  etch  or  just  an  etch.  This  includes  the  nature  and  order  in  which  chemi¬ 
cals  were  used;  and  (3)  the  amount  and  type  of  native  oxide  which  is  either 
formed  or  allowed  to  grow.  All  three  of  these  have  been  shown  to  affect  and 
alter  interface  charge  so  that  in  effect  there  is  no  generic  Si02: HgCdTe 
interface.  The  properties  of  these  interfaces  are  comparable  one  to  another 
only  when  the  surface  preparation  and  manner  of  deposition  of  the  Si02  is 
similar.  For  that  reason  and  in  order  to  form  a  meaningful  data  base  we  have 
adopted  for  this  study  a  specific  surface  process  and  deposition  procedure 
which  is  used  for  all  our  MIS  devices. 

The  interface  that  results  shows  a  very  small  voltage  hysteresis  compared 
to  the  other  charge  sources  so  that  slow  interface  charge  appears  to  be  a 
second  order  effect.  Wie  interface  trap  contribution,  as  discussed  in  the 
appendix  to  this  section,  can  be  the  dominant  source  of  interface  charge. 
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This  has  extremely  important  iaplieationa  for  the  way  in  which  SiOj  behaves  as 
a  passivation  in  diode  applications.  Specifically,  interface  charge  can  be 
expected  to  be  temperature  and  doping  level  dependent. 

Interface  trap  distributions  are  believed  to  be  due  principally  to  scat¬ 
tering  induced  tailing  of  the  conduction  and  valence  bands  into  the  gap  at  the 
surface.  1b  first  order  then  we  can  expect  these  to  be  essentially  donor  like 
(positively  charged  when  unoccupied)  and  acceptor  like  (negatively  charged 
when  occupied).  Ihese  trap  densities  are  high  at  the  band  edges  but  decrease 
by  several  orders  of  magnitude  by  mid-gap  so  that  the  individual  aid-gap 
contributions  seen  to  give  the  values  in  the  range  of  1010  to  1012  cm”2  eV1. 


that  are  measured,  this  is  illustrated  schematically  in  Figure  1 .  The  most 
important  aspect  for  determining  the  amount  of  charge  under  these  conditions 
then  is  not  so  much  the  absolute  value  of  Dit  but  its  symmetry.  A  skewed 
distribution  will  mean  that  the  higher  values  toward  one  band  edge  will  not  be 
as  effective  in  cancelling  the  opposite  charge  of  the  other.  Accurate  calcu¬ 
lation  of  this  effect  unfortunately  depend  on  details  of  Dit  in  the  region 
where  it  is  least  well  characterised.  However,  meaningful  estimates  of  the 
size  of  the  effect  can  be  made,  as  in  the  appendix,  and  do  indicate  the  amount 
of  Qit  is  in  the  appropriate  range  to  cause  observed  flat  band  shifts.  For 
use  as  a  passivation  then  an  photovoltaic  devices  the  net  charge  on  the  n  and 
p  side  will  be  different.  A  process  optimized  for  one  side  will  not  be  neces¬ 
sarily  be  optimized  for  the  other.  The  SiOj  will  act  differently  depending  on 
the  surface  fermi  energy.  This  is  illustrated  in  Figure  1  where  contributions 
to  interfaces  trap  charge  on  n-  and  p-type  surfaces  are  schematically  repre¬ 
sented.  It  becomes  clear  then  that  for  effective  passivation  it  is  important 
to  produce  very  low  and  very  symmetrical  distributions  of  interface  traps. 
Symmetrical  densities  of  less  than  1  x  101®  cm-1 ev-1  have  been  achieved  on  a 
handful  of  wafers  so  that  it  is  not  out  of  the  question.  Evidence  also  sug¬ 
gests  this  will  greatly  benefit  1/f  noise  levels  as  well  as  surface  recombina¬ 
tion  velocity. 


•  Tm  . Vj  _»  i 
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Figure  1 .  Schematic  Representation  of  Contributions 

to  Interface  Trap  Charge  on  n-  and  p-type  Surfaces 
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ELECTRICAL  PROPERTIES  OF  THE  Si02:HgCdTe  INTERFACE* 

J.A.  Wilson  and  V.A.  CoCton 
Santa  Barbara  Research  Center,  Goleta,  CA  93117 
By  determining  the  dependence  of  flatband  bias  on  insulator  thicknesses 
in  MIS  devices  made  with  PHOTOX"  Si02  on  n-type  HgQ^CdQ  2Te  we  have  for  the 
first  time  determined  the  amount  of  interface  charge  and  the  effective  work 
function  of  the  semiconductor  surface.  Flatband  bias  was  determined  from  the 
surface  potential  as  obtained  by  Berglund  integration  of  quasi-static  capaci¬ 
tance  measurements  at  77K.  Insulator  thickness  was  varied  from  400  to 
4000A.  The  amount  of  fixed  charge  was  seen  to  vary  between  -0.3  x  10^  cm-^ 
and  + 8  x  10*®  cm~^  depending  on  the  type  of  deposition  sequence  used.  It  was 
also  seen  that  H20  could  be  absorbed  into  the  layer  and  removed  by  vacuum 
baking  at  temperatures  of  60  to  100*C  and  resulted  in  a  change  of  interface 
fixed  charge  of  up  to  8  x  10*°  cm-^ .  This  is  in  agreement  with  the  range  of 
charge  variation  seen  above.  The  presence  of  large  amounts  of  H20  also  was 
seen  to  introduce  surface  states  at  0.06  and  0.14  eV  above  the  valence  band 
edge.  (Eg  »  0.26  eV  at  77K).  The  effective  work  function  difference  between 
the  Ti,  Au  bilayer  gate  and  the  HgCdTe  surface  was  determined  to  be 
♦0.03V  ±0.05V. 

INTRODUCTION 

PHOTOX"  Si02  has  become  an  important  material  for  use  in  passivating 
HgCdTe  photovoltaic  devices  in  recent  years.  The  net  charge  associated  with 
these  types  of  interfaces  is  the  dominant  parameter  determining  its  effective¬ 
ness.  PHOTOX"  is  a  form  of  low-temperature  chemical  vapor  deposited  (LTCVD) 
Si02.  These  types  of  Si02  are  achieved  by  reacting  silane  ( SiH 4)  with  oxygen 

^Supported  by  DARPA  Contract  MDA-903-83-C-0108. 
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where  different  processes  use  different  techniques  for  obtaining  the  oxygen. 
The  PHOTOZ"  process  uses  Hg  vepor  excited  by  exposure  to  UV  (2547 A)  radiation 
which  catalyzes  the  decomposition  of  nitrous  oxide  (^O)  to  produce  oxygen 
with  H2O  as  one  of  the  by-products.*  The  emphasis  in  exploring  these  process¬ 
es  is  that  the  HgCdTe  alloy  is  intolerant  of  temperatures  in  excess  of 
100*C.  Due  also  principally  to  their  low  deposition  temperatures,  the  Si(>2 
which  results  can  be  very  different  from  the  Si<>2  grown  on  Si  by  a  higher 
temperature  process. 

EXPERIMENT 

The  interface  with  HgQ.yCdg^Te  has  been  characterized  using  MIS  capaci¬ 
tors  made  on  wafers  cut  from  ingots  grown  by  the  horizontal  zone  melt  tech¬ 
nique.  These  wafers  were  cut  parallel  to  the  long  axis  of  the  ingot  without 
regard  to  crystal  orientation.  They  were  then  isothermally  annealed  in  Hg 
vapor  to  lower  the  cation  vacancies  to  less  than  5  x  10  cm  J  so  that  the 
final  electrical  activity  is  determined  by  residual  impurities.  This  is 
usually  n-type  in  the  range  of  1  *  10*^  to  1  x  10*^  cm“^.  Since  only  the 
outer  4-5  mils  are  affected  many  commonly  present  acceptors  which  are  also 
fast  diffuses  in  HgCdTe  are  gettered  to  the  unannealed  wafer  core  and  do  not 
contribute  to  the  surface  region  electrical  properties. 

Capacitors  are  then  made  by  depositing  PHOTOX”  Si02  followed  by  a  gate 
metal.  The  gate  is  done  by  thermal  evaporation  through  a  shadow  mask  and 
consists  of  a  thin  Ti  layer  (400A)  followed  by  5000A  of  Au.  Ground  contact  to 
the  n-type  HgCdTe  is  made  with  In  solder. The  insulator  thickness  was  varied 
from  500A  to  4000A  in  four  steps.  These  sets  of  devices  were  then  character¬ 
ized  by  high  (1M  Hz)  and  low  frequency  (  1C’2  Hz)  capacitance  measurements. 

The  data  was  reduced  to  give  flat  band  bias  and  interface  trap  density.  The 
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flat band  bias  is  expected  to  vary  in  a  characteristic  manner  with  insulator 
thickness  depending  on  the  distribution  of  net  interface  charge.  For  the  case 
of  charge  residing  in  a  sheet  near  the  interface  the  dependence  is 


Q*.  +  Q._ 
^fix  ^it 

C. 

ins 


where  the  insulator  capacitance  is  given  by 


e  is  the  dielectric  constant  of  the  insulator  [5  x  10  ^  —  ]  ,  A  is  the  area 
(2.04  x  10~3  cm*)  and  d  is  the  thickness.  Qf£x  and  Q^t  are  the  fixed  charge 
and  interface  trap  charge.  Taken  together  these  are  the  net  interface 
charge.  is  the  effective  work  function  difference  between  the  metal  gate 

and  the  semiconductor.  We  see  that  for  this  case  we  have  a  linear  dependence 
of  Vfj,  on  net  charge.  By  measuring  at  different  thicknesses  the  net 
<  harge  can  be  measured  from  the  slope  and  the  effective  work  function  differ¬ 
ence  from  the  zero  thickness  intercept.  Distributions  of  charge  other  than  a 
sheet  at  the  interface  will  of  course  show  nonlinear  dependence.3 

The  flatband  bias  is  determined  from  plots  of  surface  potential  (Ys) 
versus  applied  bias  obtained  by  Berglund  integration  of  quasi-static  C-V 
curves^  These  are  compared  to  measurements  of  the  shift  of  high  frequency 
curves  from  ideal  where  no  distortion  of  the  high  frequency  curve  in  depletion 
occurs . 

The  insulator  was  deposited  by  two  significantly  different  procedures. 

In  one  case  a  single  wafer  was  coated  with  successive  thin  (  500A)  depositions 
of  Si(>2  on  successively  smaller  areas  of  the  wafer.  The  first  deposition 
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covered  the  entire  wafer  end  defined  the  interface  region  for  all  different 
thickness  devices.  Thicker  devices  then  built  up  with  up  to  four  Si(>2  layers 
each  separately  deposited.  The  other  procedure  used  separate  smaller  wafers 
cut  from  a  larger  single  wafer,  each  separately  coated  with  a  single  layer  of 
appropriate  thickness.  The  devices  produced  in  this  way  are  different  in  that 
the  first  layer  and  the  interface  have  seen  process  steps  and  conditions  which 
the  single  layer  devices  have  not.  Specifically  they  were  exposed  to  the  high 
temperatures  (100*C)  and  vacutn  (1  x  10"J  torr)  of  the  predeposition  PHOTOX" 
chamber. 

RESULTS 

The  results  for  these  two  types  of  layers  show  interesting  differences. 

As  shown  in  Figure  1,  the  data  falls  into  two  groups,  both  of  which  are  reas¬ 
onable  fits  to  linear  thickness  dependence.  The  multilayer  set  of  samples 

o 

shows  a  small  positive  slope  indicating  a  net  charge  of  -3  *  107  cm  .  The 

10  —2 

single  layer  samples  show  a  larger  positive  charge  of  -*-8  x  10  cm  .  The 
zero  thickness  intercepts  for  both  sets  of  data,  are  essentially  equivalent 
considering  the  scatter  in  the  data  for  the  assumption  of  linear  thickness 
dependence.  Thus  it  is  clear  that  the  charge  at  these  interfaces  is  well 
characterized  by  a  sheet  very  near  the  actual  interface. 

The  amount  of  interface  trap  charge  contributing  to  the  net  charge  can  be 
estimated  from  the  apparent  interface  trap  density  (D^t)  derived  by  taking 
either  the  difference  between  high  and  low  frequency  C-V  curves^  or  the  dif- 
ference  between  experimental  and  ideal  curves  .  Results  are  shown  in  Figure 
2.  Also  shown  are  estimated  contributions  from  scattering  induced  tailing  of 
the  conduction  and  valence  bands  into  the  surface  energy  gap.  States  in  the 
upper  half  of  the  gap  are  taken  as  donorlike  (positively  charged  when  unoccu- 
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pied)  while  Chose  in  Che  lower  half  as  accepcorlike  (negatively  charged  when 

occupied).  An  upper  limit  Co  Che  neC  interface  Crap  charge  ac  flatband  is 

Chus  seen  to  be  on  Che  order  of  -1  *  10  cm  .  This  is  a  small  negative 

charge  on  Che  order  of  chat  seen  in  Che  multilayer  sample  and  chus  can  be  a 

significant  contribution  Co  layers  low  in  absorbed  water.  This  method  of 

estimating  interface  trap  charge  is  intended  only  to  indicate  an  order  of 

magnitude.  Further,  this  value  should  be  taken  only  as  an  upper  limit,  since 

the  existence  of  any  nonuniformities  in  either  Qfix  or  Q£t  will  make  the 

apparent  interface  trap  density  appear  artificially  high.  These  nonuniformi- 

ties  also  produce  structure  in  D£t,  especially  in  the  region  of  g .  This  is 

the  source  of  the  drop  off  in  D^t  seen  near  the  band  edges  in  Figure  2.** 

The  difference  in  deposition  sequence  between  the  two  types  of  devices  is 

significant  in  a  way  which  can  account  for  the  difference  in  net  fixed 

charge.  For  the  multilayer  samples  the  chamber  was  opened  to  move  the  shutter 

then  reevacuated  (~1  x  10-^  torr)  for  the  next  run.  This  gives  each  layer  a 

15  minute  vacuum  bake  at  100*C.  The  single  layer  samples  do  not  see  this 

step. He  have  found  that  exposing  single  layer  (1500A  thick)  MIS  devices  to  1^0 

vapor  (901  humidity,  24  hours  at  32 *C)  and  then  vacuum  bakes  (on  the  order  of 

60*C  for  12  hours)  can  shift  the  capacitance  voltage  curves  by  about  -0.8V, 

corresponding  to  a  change  in  net  charge  of  +8  x  10*®  cm~^,  equivalent  to  the 

difference  in  the  two  slopes  in  Figure  1.  This  is  illustrated  in  Figure  3. 

These  vacuum  bakes  are  much  longer  than  those  associated  with  predeposition 
* 

system  purging  but  they  are  also  at  a  significantly  lower  temperature  and  done 
on  much  thicker  layers  (1500A  compared  to  500A).  The  baked  devices  were  also 
completed,  with  metal  gates  deposited  on  the  Si02.  While  the  first  layer  of 
the  multilayer  samples  have,  at  this  stage,  an  uncoated  surface.  This  sug¬ 
gests  that  water  does  play  a  major  role  in  the  net  charge  seen  in  PHOTOX" 
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SiOj.  Although  eta*  exact  Mchtsiia  is  to  b*  determined,  water  has  been  seen 
to  play  a  maker  of  important  roles  in  PHOTO!"  ia  addition  to  the  effects  on 
net  chart*.7 

Specifically  the  presence  of  H2O  due  to  high  humidity  exposure  reversibly 
alters  the  interface  trap  distribution  (D^t)  as  illustrated  in  Figure  4.  This 
change  in  D^c  is  always  associated,  in  these  samples,  with  the  above  stent ioned 
shift  in  Vg|,  by  about  on*  volt  toward  the  negative  corresponding  to  the  posi¬ 
tive  change  ia  net  charge.  Assigning  this  change  as  due  to  a  change  in  Q^c 
requires  the  nature  of  D^t  in  the  region  of  the  peaks  at  0.06  and  0.12  eV  to 
change  frost  acceptor  to  donor  like.  These  are  of  course  filled  and  therefore 
neutral.  However  there  are  then  less  filled  acceptor  like  states  (negative 
charge)  to  compensate  the  unfilled  donor  like  states  (positive  charge)  above 
E£s  resulting  in  the  increase  in  positive  contributions  to  Q^t. 

The  intercept  with  the  vertical  axis  in  Figure  1  is  a  measure  of  the  . 
effective  work  function  difference  between  the  HgQ  >7cd()  ^  #urface  and  the 
Ti/Au  gate  we  see  that  l^8  *  +0.03V  t  0.05V.  The  surface  thus  has  very  nearly 

Q 

the  sane  work  function  as  Ti  (4.4  eV)  and  is  consistent  with  a  previous 

o 

measurement  made  at  the  anodic  oxide  interface  by  Nemirosky.  These  wafers 
are  n  type  at  1  x  10**  cm-'* ,  the  bulk  ferai  energy  is  0.213  eV  and  the  energy 
gap  at  77!  is  0.243  eV.  So  that  within  our  experimental  tolerances  the  elec¬ 
tron  affinity  for  this  surface  is  essentially  4.4  eV. 


CONCLUSIONS 

It  is  clear  that  at  the  PHOTOX"  SiOj  HgCdTe  interface  the  net  interface 
charge  (interface  trap  charge  plus  fixed  charge)  residues  in  a  sheet  very  near 
the  interface.  The  amount  of  net  charge  can  vary  through  changes  in  the 
nature  of  charged  interface  traps.  The  actual  fixed  charge  component  appears 


TAP41300 


13 


to  be  a  very  mall  amount  of  negative  charge  on  the  order  of  -1  *  10 
Finally  the  Hg^  7CdQ  3Te  electron  affinity  ia  estimated  at  4.4  eV. 
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Figure  3.  Effect  of  1^0  on  high  frequency  capacitance  measure¬ 
ment.  Shows  "water  knee"  (arrows)  and  resulting 
flatband  bias  shift. 


1.  J.W.  Paters,  N.H.  Rogers,  J.T.  Hall,  B.M.  Yes  and  D.R.  Rhiger .  Paper  No. 
202  Electro-chemical  Society  Meeting,  9-14  May  1982,  Montreal,  Canada. 

2.  J.  Tregilgas,  J.  Beck  and  B.  (Shade  "Type-conversion  of  (HgCd)Te  Induced  by 
the  Redistribution  of  Residual  Acceptor  Impurities”.  Paper  delivered  at 
the  1984  O.S.  Workshop  on  the  Physics  and  Chemistry  of  HgCdTe,  May  15-17, 
1984,  San  Diego,  CA. 

3.  A.S.  Grove.  Physics  and  Technology  of  Semiconductor  Devices,  J.  Wiley  & 
Sons,  New  York  (1967)  p.  281. 

4.  C.N.  Berglund,  tbrr  Trans.  Electron  Devices,  H3-13,  701  (1966). 

5.  R.  Cartagne  and  A.  Vapaille,  Surface  Science  28,  557  (1971). 

6.  R.  Cartagne  and  A.  Vapaille,  Surface  Science  21 ,  (197  )  157-193. 

7.  J.A.  Wilson  and  V.A.  Cotton  "Effects  of  H20  on  the  Si02/HgCdTe 

interface".  Submitted  to  J.  Appl.  Phys. 

8.  H.B.  Michaelson,  IBM  J.  Res.  Dev.,  22,  72  (1978). 

9.  Y.  Nemirovsky  and  I.  Ki-dron  Solid  State  Elect.,  22,  831  (1979). 


TAP 41 300 


i  J  m 


:vwvv. 


v.viv  v>,vw. 


Section  II 

Channeled  Implantation 


Section  II-1 

As  discussed  in  the  introduction  to  the  first  interim  report,  there  is  a 
distinct  difference  in  the  effects  on  device  performance  of  interface  charge, 
depending  on  whether  the  device  involves  ion  implantation  or  not.  These  data 
taken  on  the  effects  of  gate  bias  on  the  junction  region  of  LPE-grown  double 
layers  were  compared  to  ion  implanted  devices.  For  nominally  similar  devices 
the  RQA  of  implanted  junctions  is  lower  them  double-layer  junctions.  MIS 
devices  on  B+  implanted  surfaces  have  been  seen  to  exhibit  higher  thermal 
generation  rates  than  unimplanted  but  similarly  doped  material.  It  appears  a 
deep  state  is  introduced  by  ion  implantation  which  adversely  affects  the 
surface  region  and  therefore  device  performance.  It  is  well  documented  that 
the  resultant  n-type  doping  resulting  from  ion  implantation  into  HgCdTe  is  a 
result  of  electrically  active  deunage  sites  so  that  the  nature  of  the  specific 
ion  is  largely  irrelevant.  These  effects  occur  for  random  implants  where  each 
ion-lattice  collision  can  transfer  a  large  portion  of  the  ion  kinetic  en¬ 
ergy.  It  is  possible  to  implant  into  the  (110)  face  where  the  zinc  blend 
structure  presents  open  channels.  In  this  case  the  probability  of  direct 
collision  with  a  lattice  site  except  at  the  crystal  surface  is  low.  Ions 
channeled  in  this  way  will  lose  energy  to  the  lattice  by  interaction  of  their 
electron  clouds.  The  “nuclear"  component  of  energy  loss,  believed  responsible 
for  the  lattice  damage,  will  only  be  probable  very  near  the  end  of  the  ion 
range  where  the  ion  kinetic  energy  is  low  and  less  likely  to  cause  damage. 

This  technique  is  being  pursued  as  a  means  to  produce  electrically  active 
regions  without  the  usual  n-type  deunage.  Of  particular  interest  are  ions  of 
isotopes  of  H  and  He.  These  have  proven,  however,  very  difficult  to  detect 
once  implanted.  In  this  section  we  report  a  demonstration  of  a  technique  used 
to  detect  these  ions  in  HgCdTe  and  CdTe. 
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Section  II-2 


Depth  distributions  of  ^H,  2Hr  and  *He  implanted 
into  HgCdTe  and  CdTe  measured  by  secondary 
ion  mass  spectrometry 

R.  G .  Wilson,  Hughes  Research  Laboratories,  Malibu,  CA  90265 

Abstract 

The  capability  to  measure  atom  densities  of  ^H,  2H,  and  *He 
in  HgCdTe  and  CdTe  above  about  1018  cm“3  in  background- subtracted 
cesium  beam  secondary  ion  mass  spectrometry  profiles  has  been 
demonstrated  using  ion  implantation  and  implanted  fluences . 

Range  parameters  are  determined  for  ion  energies  between  150  and 
300  keV.  Projected  ranges  for  %  are  slightly  greater  than  those 
for  2H  in  both  HgCdTe  and  CdTe.  Hydrogen  ion  ranges  are  about 
0.75  ym/100  keV  of  energy  in  both  materials,  and  the  third  moment 
of  the  depth  distributions  is  about  -2.1. 
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In  1971,  Foyt ,  Harman,  and  Donnelly1  reported  the 
fabrication  of  n-p  junction  photodiodes  in  Hg1-J{cdxTe  (x  «  0.2- 
0.5)  by  the  implantation  of  protons  of  100  keV  energy  and  1  to  5 
x  1014  cm"2  f luences .  Electrical  performance  parameters  were 
reported,  which  included  the  observation  that  the  n  layer  created 
by  the  proton  bombardment  was  about  1  urn  thick,  but  no 
information  about  the  depth  of  penetration  of  the  100-keV  protons 
or  the  depth  profile  or  density  of  the  implanted  hydrogen  atoms 
was  given.  No  work  has  been  subsequently  reported  about  the 
implantation  of  hydrogen  or  helium  ions  in  HgCdTe  although 
considerable  work  has  been  reported2"14  using  heavier  ions, 
especially  boron.2"11  Of  the  10  papers  that  used  B  implantation, 
four  measured  the  B  depth  distribution  by  using  secondary  ion 
mass  spectrometry  (SIMS)  profiling.2"5  N-type  regions  are 
frequently  reported  to  be  produced  by  implantation  (damage) 
irrespective  of  which  ion  is  implanted,  and  the  n  regions  are 
deeper  than  the  implanted  ion  ranges . 

In  1967,  Smirnov  and  Kachurin  et  al.15»16  reported 
electrical  effects  of  implanting  Cd  and  In  intoCdTe.  Donnelly 
et  al  *12  reported  t:ype  conversion  and  the  formation  of  p-n 
junctions  by  the  implantation  of  400-keV  As  into  nrCdTe. 
Subsequently,  others  reported19-20  the  results  of  implanting 
heavy  ions  into  CdTe.  No  work  has  been  reported  in  which  the 
light  ions  H  or  Be  have  been  Implanted  into  CdTe. 

Implantation  ranges  are  shallow  for  heavy  ions  in  the  high 
atomic  number  substrates  HgCdTe  and  CdTe.  The  surfaces  of  these 
materials,  especially  HgCdTe,  are  subject  to  effects  caused  by 
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defects  and  vacancies .  Light  ions  (e  .g H  and  He)  cause  fewer 
defects  and  penetrate  more  deeply  than  heavy  ions .  Light  ions 
therefore  may  be  useful  for  creating  deeper  n-p  junctions  in 
either  HgCdTe  or  CdTe,  as  Foyt  et  al .  reported  for  100  keV 
protons  .  We  have  implanted  the  light  ions  ^H,  2H,  and  *He  into 
HgCdTe  and  CdTe  at  energies  from  ISO  to  300  keV  and  profiled 
their  depth  distributions  using  SIMS.  We  have  performed  Pearson 
IV  fits  on  these  depth  distributions  and  report  the  peak  ranges 
Rn,  the  projected  ranges  Rp  (y) ,  the  range  straggles  Ar^  (o) ,  the 
skewness  or  third  moment,  and  the  kurtosis  or  fourth  moment,  that 
is,  the  four  moments  of  the  depth  distributions,  and  the  doping 
densities  achieved.  This  is  the  first  time  that  profiles  for  H 
or  He  impurities  in  HgCdTe  or  CdTe  have  been  reported.  These 
results  also  indicate  that  hydrogen  (as  *H  or  hi)  can  be  profiled 
as  an  impurity  in  HgCdTe  or  CdTe  by  secondary  ion  mass 
spectrometry,  if  the  H  atom  density  is  greater  than  about  5  x 
101®  cm”3 . 

Implantation  into  HgCdTe  and  CdTe  was  performed  as  it  is  for 
Si,  GaAs,  and  other  semiconductors.  Samples  of  (111)  horizontal 
zone  refined  and  (111)  LPE  HgCdTe  on  CdTe  substrates  and  of  the 
semi-insulating  CdTe  substrates  themselves  were  implanted  with 
1H,  ha,  and  *He  ions  at  energies  between  ISO  and  300  keV, 
fluences  of  10**,  10*5,  and  10*6  cm*2,  and  at  current  densities 
of  <v  10~  A  an”2  at  room  temperature.  The  samples  were  tilted  7° 
from  the  surface  normal  to  minimize  channeling.  Implantation 
energies  are  accurate  to  within  about  2%  and  fluences,  to  within 
about  10%. 
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Establishing  a  SIMS  technology  for  profiling  hydrogen  and 
helium  in  HgCdTe  and  CdTe  required  special  effort  because  their 
ionization  potentials  are  high,  their  electron  affinities  are 
low,  and  H  is  present  in  the  vacuum  ambient.  Profiles  could  not 
be  obtained  for  the  101*  cm*2  implant  fluences  in  either 
material,  and  only  for  the  10*®  cm-2  He  fluence  in  HgCdTe.  We 
were  not  successful  with  He  in  CdTe  even  for  the  10*®  cm"2 
fluence,  probably  because  of  the  semi-insulating  nature  of  the 
CdTe  and  the  associated  slight  charging  difficulties.  The 
semi- insulating  CdTe  samples  were  coated  with  about  40  run  of  Au 
to  minimize  charging  during  SIMS  profiling.  *H  and  2H  were 
profiled  using  Cs  ion  bombardment  and  negative  ion  spectrometry. 
Helium  was  profiled  using  oxygen  ion  bombarduent  and  positive  ion 
spectrometry.  Sputtering  rates  of  about  10  nm/s  were  used. 
Background-subtracted  detection  limits  (atom  densities)  of  10l8 
cm”2  were  achieved  for  *H  and  *He,  and  *  10*^  cm-2  for  2H.  The 
SIMS  measurements  were  carried  out  using  a  CAMECA  IMS  3f 
instrument  at  Charles  Evans  and  Associates .  The  atom  densities 
were  determined  by  equating  the  integral  of  the  profile  to  the 
implanted  fluence .  The  depth  scales  were  obtained  by  measuring 
the  SIMS  crater  depths  using  a  surface  profilometer  (Dektak)  and 
the  crater  depths  are  accurate  to  about  +  7% . 

The  depth  distributions  of  *H  and  2H  in  HgCdTe  are  shown  in 
Figure  1  and  in  CdTe,  in  Figure  2.  The  profile  for  He  in  HgCdTe 
is  shown  in  Figure  3.  Pearson  IV  fitting  was  performed  on  all 
profiles.  The  values  of  the  maximum  atom  density  Nmax,  depth  at 
the  maximum  R^,  the  projected  range  Rp  (or  u),  the  range  straggle 


ARp  (or  <T)  ,  the  skewness  Y^,  the  kurtosis,  B2,  and  the  relative 
kurtosis  82  are  given  in  Table  I.  The  profiles  for  and 
were  measured  in  the  same  SIMS  profile,  so  their  relative  depths 
are  subject  to  the  2%  energy  inaccuracy,  not  the  +  7%  crater 
depth  inaccuracy.  The  range  of  or  2H  in  either  material  is 
approximately  0.75  urn  per  100  keV  of  H  ion  energy.  The  third 
moments  for  H  in  both  materials  are  moderately  negative  -2) , 
consistent  with  theory  for  light  ions  in  high  atomic  number 
substrates.  The  statistics  were  poorer  in  the  data  for  H  in 
CdTe,  so  we  show  the  actual  data  points  and  how  we  chose  to  draw 
the  depth  distribution  curve  through  them  in  Figure  3.  The 
statistics  for  H  in  HgCdTe  were  good  enough  to  clearly  define  the 
curves  shown  in  Figure  1. 

Calculations  of  projected  range  Rp  have  been  made  for  H  in 

CdTe  by  Gibbons  et  al  .2*  using  the  Lindhard,  Scharff,  and  Schi$tt 

formalism.  Their  calculated  values  of  R  and  ARn  for  150  keV  H 

P  P 

in  CdTe  are  1.40  and  0.30  ym,  respectively;  our  experimental 
values  for  HgCdTe  are  1.11  and  0.31  ym,  respectively.  Their 
values  for  300  keV  H  in  CdTe  are  2.78  and  0.40  urn,  respectively; 
our  experimental  values  for  CdTe  are  2.15  and  0.33  urn  for  ^H  and 
1.96  and  0.44  ym  for  2H,  and  for  HgCdTe,  our  experimented  values 
are  2.35  and  0.64  ym  for  ^H.  The  differences  between  calculated 
and  experimental  values  is  within  20  or  30%.  This  is  very  good 
agreement  for  the  large  difference  in  atomic  number  between  ion 
(1)  and  substrate  (y  50)  . 

Foyt  et  al  reported  an  n  region  created  in  HgCdTe  by  100 
keV  proton  implantation  that  was  about  1  ym  deep.  He  have  shown 


here  that  £or  100  keV  protons  the  Rp  of  the  distribution  should 
be  about  0.75  ym  and  that  the  depth  distribution  would  be  about  1 
m  deep  about  one  decade  below  the  peak  of  the  depth  distribution 
on  the  deep  side  of  the  distribution,  in  good  agreement  with  the 
n  region  reported  by  Foyt  et  al .,  assuming  that  carrier  removal 
and/or  doping  effects  of  defects/ dam age  are  caused  by  implanted 
protons  at  densities  down  about  one  decade  from  the  value  of  the 
peak  of  the  implanted  proton  depth  distribution,  which  should  be 
between  2  x  10^a  and  1  x  10^  cm"  3  based  on  their  fluence  of  1  to 
5  x  1014  cm" ^  and  our  experimental  profiles. 

The  work  was  supported  by  the  Materials  Sciences  Office  of 

^{IT  Aiv y^-C-OI  0?. 
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Figure  Captions 


Figure  1. 

Figure  2. 


Depth  distributions  (measured  using  SIMS)  o£ 
hydrogen  implanted  into  HgCdTe . 

Depth  distribution  (measured  using  SIMS)  of  helium 
implanted  into  HgCdTe . 

Depth  distributions  (measured  using  SIMS)  of 
hydrogen  implanted  into  CdTe. 
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Section  III 
Surface  Study 


The  surface  and  interface  properties  of  Hg1 _xCdxTe  cannot  be  understood 
without  consideration  of  such  conventionally  bulk  properties  as  the  nature  of 
the  bonds  which  make  up  the  lattice  or  the  types  and  motion  of  extended  de¬ 
fects  in  the  material.  The  width  of  the  liquid  phase  epitaxially  grown  (LPE) 
HgCdTe/CdTe  he tero junction  and  accompanying  broad  distribution  of  disloca¬ 
tions1  are  examples  of  the  way  in  which  driving  forces  such  as  the  difference 
in  cation  bonding  can  affect  what  is  normally  considered  interfacial  beha¬ 
vior.  The  relationship  between  out-diffusion  and  evaporation  of  Hg  from  the 
surface  to  sample  history  or  surface  preparation  similarly  indicates  the  need 
for  a  broadly  based  understanding  in  order  to  fully  comprehend  and  solve 
surface  and  interface  related  problems. 

This  outlook  is  reflected  in  the  work  reported  here.  The  electronic 
structure  of  the  alloy  at  once  accounts  for  the  desired  tunable  bandgap  of  the 
materials  as  well  as  many  materials  problems  in  growth  and  preparation2.  A 
detailed  description  of  the  electronic  structure  can  be  gained  through  the 
technique  of  angle-resolved  photemission  spectroscopy;  a  preliminary  report  on 
this  project  is  the  subject  of  Section  III-2.  Section  III-3  presents  results 
obtained  using  transmission  electron  microscopy  to  characterize  extended 
defects  observed  at  the  LPE  HgCdTe/CdTe  he tero junction.  This  information  is 
of  value  not  only  in  characterizing  an  interface  of  extreme  technological 
importance  but  also  because  analysis  of  the  imaged  defects  can  yield  data  on 
stacking  fault  energies  and  dislocation  motion  and  annealing  -  data  which 
strengthens  the  link  between  the  bonding  description  originating  in  the  elec¬ 
tronic  structure  studies  and  macroscope  materials  properties. 

To  supplement . the  perspective  gained  from  studies  of  the  bonding,  the 
surface  properties  were  probed  directly  in  a  comparative  study  of  activated 
oxygen  uptake  on  the  clean  (110)  faces  of  CdTe,  HgTe,  and  HgQ  ggCdg  31^  an<* 
in  a  first  attempt  at  learning  more  about  the  mechanism  of  surface  type  con¬ 
version  observed  upon  cleavage^  by  performing  the  cleaving  at  reduced  tempera¬ 
ture.  These  efforts  are  discussed  in  Sections  III-4  and  III-5  respectively. 

By  an  appropriate  modification  to  the  experimental  procedure  discussed  in 
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Section  III-4,  initial  data  regarding  the  kinetics  of  the  oxidation  process 
were  obtained  which  augment  information  on  the  chemistry  of  the  oxygen-HgCdTe 
interaction.  As  pointed  out  in  this  section,  both  the  chemical  and  kinetic 
trends  cam  be  understood  in  the  framework  of  the  alloy  bonding. 
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Section  III-2 


II.  Angle-Resolred  Photo  ealseion  Spectroscopy  of  Hg^Cd^Te: 
Preliminary  Analysis  of  data  froa  CdTo. 

A)  Introduction. 

Ones  tho  oloetronie  structure  of  a  sealconductor  lo  accurately  known, 
proportion  ranging  froa  tho  conduction  hand  oloctron  effective  aaoo  to  tho 
•tacking  fault  energy  can  bo  ontiaatod.  In  tho  cano  of  Mercury  Cadaiua 
Tollurido  CUCT) ,  tho  rolationnhip  botvoon  aatorialo  proportieo  and 
oloetronie  otructuro  in  particularly  evident .  Bocauno  of  tho  higher 
nuclear  charge  of  Hg  coapared  to  Gd.  tho  seloctire  penetration  of  tho  e- 
oloctron  wave  function  to  the  aucloun.  and  tho  influence  of  tho  junt 
conplotod  d10  shell,  tho  binding  energy  of  tho  atoaic  Hg  6n  valence 
electrons  in  1.4  oV  groator  than  that  of  tho  Cd  6n  levels  while  tho 
(unoccupied)  ralonco  p-orbitaln  hare  niailar  onorgion.  Thin  difforonco  in 
tho  atoaic  ralonco  Staten  is  central  to  tho  alloy  oloetronie  otructuro  and 
is  strikingly  aanlfost  in  two  aspects  of  tho  solid  ralonco  band  otructuro 
when  Hg  and  Cd  are  distributed  in  tho  cation  sublattico  of  tho  MCT  alloy.1 

On  tho  one  hand,  tho  states  which  coaprlso  tho  conduction  band  ainiana 
arise  froa  tho  cation  ralonco  s-electrons;  as  Hg  is  substituted  for  Cd. 
tho  rolatiristic  interactions  of  tho  high  Z  Hg  core  which  are  responsible 
for  tho  increase  in  atoaic  ••oloctron  ionisation  potential  result  in  tho 
decrease  in  fundaaontal  bandgap  of  tho  alloy  with  increasing  Hg  content.2 
Thus  tho  aatorialo  property  which  aakos  MCT  attractiro  in  practical 
applications,  tho  ability  to  tune  tho  optical  absorption  edge,  is  cloarly 
a  result  of  tho  unique  electronic  structure  of  tho  alloy.1 

A  second  consequence  of  tho  atoaic  ralonco  structure  of  Hg  and  Cd  for 
tho  alloy  electronic  states  relates  to  tho  bonding  of  Hg  and  Cd  in  tho 
lattice.1  In  other  soaiconductor  alloys  studiod  to  date,*'4  tho  rlrtual 
crystal  approxiaation  (VGA)  proridos  an  adequate  description  of  the  alloy 
ralonco  bands.  In  tho  VCA.  use  is  Bade  of  a  coapositionally  weighted 
arorago  of  tho  potentials  of  tho  distributed  species  —  A1  and  Ga  in 
Alj.gOa^As.  for  oxaaplo.  This  procooduro  allows  ono  to  treat  tho  randoa 
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alloy  with  technique*  applicable  to  a  crystal  with  a  periodic  potential; 
however,  it  a  priori  renders  the  bonds  between  the  various  cations  and 
anions  indistinguishible.  In  the  ease  of  MCT.  the  difference  in  potential 
between  Hg  and  Cd  is  sufficiently  great  that  a  description  based  the  VCA 
fails  for  the  cation  s-electron  derived  states  deep  in  the  valence  bands 
3.5  to  8  eV  below  the  valence  band  maximum.5  It  is  clear  that  as  a  result 
of  the  distinct  Hg  or  Cd  contribution  to  the  alloy  electronic  structure, 
the  Hg-Te  and  Cd-Te  bonds  in  MCT  are  different.  The  consequences  of  this 
bonding  difference  for  properties  such  as  dislocation  formation  and 
hetero junction  growth1  thus  follow  froa  the  nature  of  the  alloy  valence 
electronic  structure. 

It  should  be  noted  that  the  breakdown  of  the  VCA  is  selective  in  MCT.5  As 
mentioned  above,  it  occurs  most  strongly  for  the  s-derived  states  deep  in 
the  valence  bands.  These  states  have  k-values  near  the  Brillouin  zone 
boundary.  At  the  zone  center  and  the  valence  and  conduction  band  edges, 
the  VCA  appears  to  be  successful.  The  group  of  valence  states  with  energy 
within  about  3.5  eV  of  the  valence  band  maximum  (VBM)  similarly  exhibit 
behavior  more  nearly  that  expected  in  a  periodic  crystal  or  one  in  which 
the  disorder  in  potential  is  small.  (The  states  within  3.5  eV  of  the  VBM 
have  mainly  Te  p-electron  character  and  are  therefore  less  sensitive  to 
the  relativistic  interactions  which  characterize  the  difference  in  Hg  and 
Cd  potintial.  )  In  order  to  fully  characterize  the  selective  breakdown  of 
the  VCA  and  thereby  optimize  the  accuracy  of  calculations  of  the 
electronic  structure  and  the  materials  properties  that  follow  from 
knowledge  of  the  electron  states,  a  technique  sensitive  to  both  the 
electron  energy  and  k  vector  in  the  solid  is  required.  This  sensitivity  is 
provided  by  angle-resolved  photoemission  spectroscopy7,8  (ARPES)  in  which 
both  the  photoemitted  electron  energy  and  angle  of  emission  from  the  solid 
are  determined. 

B)  Electron  State  selectivity  of  ARPES. 

The  electron  states  in  a  solid  are  characterized  by  the  quantum  numbers 
energy,  E,  and  wave  vector,  k.  In  an  ARPES  experiment,  electron  emission 
induced  by  absorption  of  ultraviolet  light  is  collected  in  a  small  solid 
angle  and  in  a  specific  direction  relative  to  the  crystal  axes  and  energy 


analyzed.  The  energy  and  ware  rector  of  the.  free  electron  are  thne 
determined.  (The  magnitude  of  the  ware  rector  is  fixed  by  the  electron 
energy,  and  the  direction  is  determined  by  the  angular  position  of  the 
detector.)  The  ability  to  relate  these  ralnee  outside  the  crystal  and  E 
and  k  for  the  states  of  the  solid  prorides  the  utility  of  ARPES  as  a  probe 
of  the  band  structure  of  solids. 

Photoemission  can  be  riewed  as  a  three  step  process  (Figure  1) ,9 
Absorption  of  a  photon  produces  an  optical  transition  between  two  allowed 
states  of  the  solid  separated  by  the  photon  energy.  If  the  excited 
electron  has  a  reloci ty  component  in  the  direction  of  the  surface,  it  may 
propagate  toward  the  surface;  the  occurance  of  inelastic  scattering  during 
this  step  will  restore  electrons  from  the  primary  distribution.  If  the 
final  state  energy  to  which  the  electron  was  excited  is  sufficiently 
great,  the  electron  has  some  probability  of  escaping  the  crystal  once  it 
reaches  the  surface.  In  the  racuum.  the  photoelectron's  kinetic  energy 
measured  by  the  detector  is  referenced  to  the  racuum  lerel.  The  ralence 
band  maximum  surras  as  a  conrenient  zero  of  energy  in  the  crystal;  hence, 
the  measured  kinetic  energy  and  the  final  state  energy  in  the 
photoabsorption  erent  deffer  by  the  sum  of  the  semiconductor  bandgap  and 
electron  affinity  (see  Figure  1) .  Because  energy  is  conserred  in  the 
optical  excitation  process,  the  initial  state  lies  h?  below  the  final 
state  in  energy.  Thus,  the  measured  photoelectron  kinetic  energy  can  be 
used  to  determine  respectirely  the  ralues  E;  and  EA  of  the  conduction  band 
and  ralence  band  states  inrolred. 

When  combined  with  conserration  of  energy,  conserration  of  momentum  in 
the  optical  excitation  stop  of  the  photoemission  process  greatly  restricts 
the  number  of  possible  transitions  between  ralence  and  conduction  band 
states  in  the  solid.  Figure  2  illustrates  a  reciprocal  space 
representation  of  optical  transitions  between  hypothetical  initial  and 
final  state  bands  for  two  photon  energies.  In  a  categorization  of  the 
electron  state  ware  rectors  in  which  all  k-ralues  are  reduced  to  the  first 
Brillouin  zone,  conserration  of  momentum  implies  the  initial  and  final 
states  must  hare  the  same  ralue  of  k.  While  the  transition  indicated  by 
the  solid  arrow  at  kt  which  connects  states  of  the  same  k  ralue  separated 
in  energy  by  h^  may  occur,  the  transition  to  the  same  final  state  from 
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the  initial  state  at  k2  (dashed  arrow  labeled  ht^)  satisfies  energy 
conservation  but  not  conservation  of  momentum  and  thus  is  not  strictly 
allowed.  The  transition  illustrated  by  the  vertical  arrow  is  refered  to  as 
a  direct  transition  whils  the  absolption  from  k  equal  k2  is  nondirect. 
Nondirect  transitions  do,  in  fact,  occur8,  as  will  be  evident  in  the  data. 
As  the  photon  energy  is  varied  the  optical  excitation  will  occur  from 
states  at  different  Et  and  k.  (See  hv2  in  Figure  2.)  The  variation  in  Et 
and  Ef  with  h v  will  result  in  a  shift  in  the  kinetic  energy  of  the 
photoelectron  peak  in  the  measured  spectrum  characteristic  of  the 
dispersion  of  the  valence  and  conduction  bands  involved.  In  addition  to 
conservation  of  energy  and  momentum,  another  constraint  is  imposed  in  the 
photoexcitation  process.  Selection  rules  determine  for  a  given 
polarization  which  bands  may  participate  in  an  optical  transition  based  on 
the  symmetry  of  the  electron  states  which  comprise  the  band  and  the 
direction  of  the  optical  electric  field.10  Vhile  k  equal  k3  or  k4  is 
allowed  for  h^2  in  Figure  2.  for  example,  the  optical  selection  rules  may 
exclude  photoexcitation  at  k4. 

As  discussed  above,  conservation  of  momentum  in  the  excitation  step 
ensures  the  equality  of  kA  and  kf.  When  the  excited  electron  escapes  over 
the  surface  potential  barrier,  a  refraction  phenomenon  occurs  such  that 
only  the  component  of  kf  (and  ki)  parallel  to  the  surface  plane  is 
conserved.7'8  Because  the  parallel  component  of  the  free  electron  wave 
vector  outside  the  solid  is  determined  by  the  angular  position  of  the 
detector  for  a  given  electron  energy,  the  choice  of  angle  of  detection 
similarly  determines  the  value  of  k  to  which  the  experiment  is  sensitive 
inside  the  solid.  The  indeterminacy  in  the  component  of  kf  perpendicular 
to  the  surface  means  that  the  valence  and  conduction  band  states  which 
contribute  to  the  ARPES  spectrum  lie  along  a  line  in  reciprocal  space 
given  by  k  equal  to  a  constant  which  depends  on  the  electron  energy.  It  is 
the  ability  to  select  the  transitions  occuring  in  the  small  volume  of  k- 
space  along  this  line  out  of  all  possible  allowed  transitons  that  makes 
ARPES  such  a  sensitive  probe  of  electronic  structure. 

C)  Experimental 

The  ambiguity  in  determining  the  value  of  k  when  only  k||is  known  can  be 


overcome  by  cartful  analysis  of  the  photoemission  data,  aa  dit catted  below 
(••ctios  0) .  Such  analytic  it  simplif itd  for  particular  choices  of 
k.  Petitioning  of  tbo  detector  along  the  normal  to  the  eaaple  turf ace 
enauret  that  the  atatea  contributing  to  the  measured  energy  distribution 
lie  in  the  direction  of  k-space  normal  to  the  plane  of  the  surface.  In  the 
data  reported,  the  (110)  face  of  Hg^Cd^Te  (**10,  0.31.  0.23,  and  0.0) 
was  prepared  by  cleaving  in  vacuum,  and  data  vae  collected  in  the  normal 
emietion  geometry.  Thus  the  measurements  are  sensitive  to  valence  and 
conduction  band  states  in  the  (110)  direction  or  F27KSX  symmetry  axis  of 
the  Brillouin  cone  (Figure  3) .  Additional  measurements  were  made  with  kA 
chosen  to  lie  on  the  hexagonal  face  of  the  Brillouin  cone  (VLV  in  Figure 
3)  for  initial  states  of  high  s-character  4.6  to  6  eV  below  the  VBM  in 
order  to  explore  the  breakdown  of  the  VCA  for  states  at  the  cone  boundary. 
This  data  will  appear  in  a  future  report.  Spectra  were  recorded  for  photon 
energies  between  13  and  27  eV  using  the  8  degree  beam  line  at  the  Stanford 
Synchrotron  Radiation  Laboratory.  By  varying  the  angle  of  incidence  of  the 
photon  beam  between  20  and  64  degrees,  the  polarisation  could  be  similarly 
altered.  The  component  of  the  electric  field  vector  in  the  plane  of  the 
surface  was  parallel  to  the  (110)  direction. 

D)  Results  and  Method  of  Analysis  for  CdTe. 

Analysis  of  the  spectra  from  HgTe  and  the  alloy  samples  is  continuing. 
Vhile  the  analysis  of  the  CdTe  data  is  not  yet  in  its  final  form, 
sufficient  progress  has  been  made  to  demonstrate  the  ARPES  technique  for 
band  structure  determination.  Normal  emission  data  as  a  function  of  photon 
energy  from  (110)  CdTe  for  20  degree  (stronger  s  polarization)  and  64 
(stronger  p  polarization)  angles  of  incidence  are  presented  in  Figures  4 
and  6  respectively.  Both  sets  of  spectra  exhibit  similar  structure  which 
generally  disperses  as  the  photon  energy  is  varied,  indication  origin  in 
bulk  rather  than  surface  states.  Increasing  the  p  polarization  content  of 
the  incident  beam  enhances  the  structures  which  for  hi«18  eV  occur  near 
>4.6.  -3.6.  and  -1.2  eV,  eometimee  obscuring  weaker  peaks.  The  increase 
in  strength  of  these  peaks  indicates  a  character.10 

Vhile  the  use  of  normal  emission  restricts  the  region  of  k-space 
contributing  to  the  photoenission  to  the  £  axis,  the  exact  locations  along 
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this  list  where  tilt  transitions  occur  vast  bt  determined  bj  analysis.  The 
criteria  to  be  met  bj  the  proper  assignment  of  the  structure  is  the 
spectra  to  parts  of  the  Brillouin  sone  include  the  correct  measnred  final 
state  energy,  Ef.  and  the  correct  spearation  between  and  Ei  such  that 
final  and  initial  states  lie  on  the  conduction  and  valence  bands  at  the 
same  value  of  k.  One  approach11  succesfully  employed  in  other  cases 
assumes  a  free  electron  parabola  for  the  final  state  with  energy  origin  VQ 
below  the  VBM;  from  the  relation 

Ef  =  W/2m  -  V0 

one  caui  determine  k  and  Ef  and  thereby  map  the  valence  bands.  Vhen  an 
accurate  initial  state  band  structure  is  known,  a  reverse  proceedure  can 
be  used  to  determine  the  high  lying  conduction  band  final  states. 12,13  In 
this  case  Ef  is  determined  from  the  photoemission  spectra  by  first 
locating  the  position  of  the  VBM  in  the  data.  From  the  band  structure,  one 
can  then  find  the  value  of  k  (there  may.  in  fact,  be  more  than  one) 
corresponding  to  Ef  for  a  given  peak  in  the  spectrum  and  photon  energy  and 
then  locate  Ef  at  this  value  of  k. 

Band  structure  calculations  for  the  initial  statees  in  CdTe  in  the  (110) 
direction  have  been  published.14'15  These  calculations  do  not  include 
conduction  bands  with  high  enough  energy  for  direct  comparison  to  the 
data.  The  analytical  approach  adopted  was  to  assume  the  nonlocal 
pseudopotential  calculation  of  Chelikowsky  and  Cohen  (reference  14)  to  be 
correct,  determine  all  possible  kf  corresponding  to  a  photoemission  peak 
at  Ef.  and  hence  locate  candidate  final  states  at  each  kf.  A  choice  of  a 
particular  k  from  the  set  of  possible  wave  vectors  was  then  made  by 
comparing  the  candidate  Ef(kf)  with  the  set  of  free  electron  parabolas  of 
the  face  centered  cubic  lattice. 

The  selected  final  states  and  free  electron  parabolas  are  shown  in 
Figure  6.  In  most  cases,  the  final  states  are  found  to  lie  near  the 
segment  of  free  electron  parabola  which  was  translated  into  the  first 
Brillouin  zone  by  reciprocal  lattice  vector  (U0)4r/a  and  is  thus  part  of 
the  parabola  centered  at  r  .  It  should  thus  be  the  major  contributor  to  the 
photoemission  process,11  as  was  indeed  found.  The  group  of  final  states 
rising  in  energy  near  k  equal  0.85  are  due  to  transitions  arising  from  the 
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high  density  of  states  region  with  initial  atnta  energy  of  approximately 
-4.6  eV.  While  the  intenaity  of  eaieaion  in  thia  peak  ia  atronger  vhen  a 
final  atate  band  ia  available,  the  tranaitiona  nonethaleaa  occur  for  all 
photon  energiea.  Thia  ia  an  axaapla  of  a  nondirect  transition;8  provided 
the  initial  atatea  lie  along  E.  the  abaenee  of  other  bande  at  thia  binding 
energy  allova  assignment  of  the  structure  to  the  lovest  valence  band.  (See 
Figure  7) . 

To  obtain  agreement  between  the  measured  final  states  and  the  free 
electron  bands,  a  value  of  inner  potential.  VQ.  of  4.3  eV  vaa  required. 
Thia  choice  of  VQ  is  consistent  vith  that  obtained  from  analysis  of  the 
Lov  Energy  Electron  Diffraction  (LEED)  intenaity  profiles  by  Duke  el  a/18, 
who  determined  a  value  of  4  eV.  (In  GaAs.  the  ARPES  value  of  9.32  eV  ia  to 
be  compared  to  8  eV  from  LEED  data.11)  The  initial  states  contributing  to 
the  spectrum  based  on  thia  preliminary  analysis  are  shown  in  Figure  7. 

E)  Future  Analysis 

To  examine  the  accuracy  of  band  structure  calculations,  a  better 
representation  of  the  final  state  bands  is  needed.  The  reasonable 
agreement  with  the  free  electron  bands  indicates  that  it  may  be  possible 
to  calculate  the  true  conduction  band  structure  vith  merely  small 
perturbations  to  the  free  electron  parabolas.  Analysis  of  the  data  for  the 
alloys  must  also  be  performed.  The  similarity  in  the  emission  from  these 
compounds  should  aid  in  establishing  a  consistent  assignment  of  the  peaks. 
Little  use  of  the  polarization  dependence  of  the  spectra  has  been  made  so 
far.  and  this  information  will  likewise  give  confidence  to  the  electronic 
structure  determination.  Finally,  the  aspect  of  the  electronic  structure 
of  MCT  key  to  understanding  the  alloy  bonding  properties  is  the  selective 
breakdown  of  the  VCA;  examination  of  the  composition  dependence  of  the 
width  of  emission  from  the  s-electron  derived  valence  states  near  the  zone 
boundary  should  provide  information  on  parameters  unique  to  non-VCA 
calculations17,18  of  the  alloy  electronic  structure. 
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Figure  Caption* 


Figure  1. 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 
Figure  6. 


Schematic  description  of  photoemission  as  a  three  step 
process  involving  photoexcitation,  transport  of  the 
excited  electros  distribution  to  the  surface,  and  escape 
into  the  vacuum. 

Reciprocal  space  representation  of  photoexcitation  step 
illustrating  constraints  on  k  and  E  imposed  by 
conservation  of  energy  and  momentum.  (See  text.) 

First  Brillouin  sons  for  a  face  centered  cubic  crystal. 
The  lines  probed  in  the  reported  ARPES  experiment  are 
shown  dashed. 

Normal  emission  energy  distribution  curves  from  (110)  CdTe 
plotted  against  initial  state  energy  (relative  to  the 
valence  band  maximum) .  The  spectra  were  recorded  for 
13<hx^(26  eV  in  0.6  eV  steps  for  polarization  20  degrees 
from  the  surface  and  along  the  (110)  direction. 

Same  as  Figure  4  but  with  the  polarization  vector  64 
degrees  from  the  surface  along  the  (110)  direction. 

Free  electron  parabolas  for  the  body  centered  cubic 
Ttcciprocal  lattice  along  £.  The  solid  curves  are  the 
parts  of  the  free  electron  final  state  parabola  centered 
at  the  origin  and  should  be  the  strongest  contributor  to 
the  photoemission.  The  points  represent  measured  final 
states  if  the  initial  bands  of  Chelikovsky  and  Cohen14  are 
assumed. 

Initial  state  bands  of  Ref.  14  and  location  of  structure 
observed  in  the  ARPES.  The  structure  was  assumed  to  match 
the  band  calculation  in  order  to  map  the  final  states. 
(See  text.) 


Figure  7. 
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Section  III-3 

TEM  STUDIES  OF  DEFECTS 

Detailed  defect  structural  characterization  of  MCT  has  been  established 
as  an  essential  prerequisite  to  understanding  and  improving  this  material. 

The  p-  to  n-type  conversion  upon  mechanical  damage  [1-3]  or  ion  implantation, 
width  of  the  LPE  grown  MCT  on  CdTe  heterojunction  [4],  and  Hg  out-diffusion 
after  mechanical  damage  [5]  are  a  few  of  the  observed  phenomena  in  MCT  that 
link  the  mechanical  defect  structure  with  the  electronic  structure.  Trans¬ 
mission  Electron  Microscopy  (TEM)  is  probably  one  of  the  most  useful  techniques 
in  viewing  directly  the  defect  microstructure  associated  with  these  phenomena, 
and  has  been  employed  here  with  the  primary  focus  on  the  LPE  grown  heterojunc-r 
tion  and  the  mechanical  damage  resulting  from  ion  milling.  The  damage  induced 
by  ion  milling  is  believed  to  be  similar  to  that  introduced  by  ion  implantation. 

In  order  to  gain  insight  into  the  defect  structure  associated  with  the 
liquid  phase  epitaxy  (LPE)  grown  MCT  on  CdTe  interface,  TEM  samples  have  been 
prepared  for  view  in  cross-section.  Because  it  is  of  great  interest  to  view 
the  defect  distribution  as  a  function  of  distance  from  the  metallurgical  junction, 
cross-sectional  TEM  samples  were  prepared  since  it  is  possible  to  view  this 
profile  directly  in  only  one  successfully  prepared  sample.  The  heterojunction 
studied  was  a  LPE  grown  MCT  (x  *  0.3)  on  a  (111)  CdTe  substrate  supplied  by 
SBRC.  This  LPE  material  was  grown  from  a  Te  melt  to  a  layer  thickness  of 
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15.8  u  m. 


The  samples  were  prepared  by  cutting  two  adjacent  1.0  x  1.5  mm  rectangular 
pieces  from  the  LPE  substrate  and  cementing  them  together  so  that  the  LPE  MCT 
layers  were  sandwiched  between  the  CdTe  substrates.  Care  was  taken  to  insure 
that  the  (110)  plane  was  oriented  parallel  to  one  of  the  sample  edges,  so  that 
lattice  images  could  be  taken  if  the  preparation  of  the  sample  was  successful. 
This  resulting  "sandwich"  was  then  imbedded  in  Bakelight,  which  acts  as  a 
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supporting  material  when  the  sample  is  very  thin.  The  samples  were 
then  mechanically  thinned  to  a  thickness  of  approximately  150  inn,  and  then 
dimpled  so  that  approximately  10  w  m  of  material  was  left  in  the  center 
cross-section  of  the  sample.  The  final  thinning  step  employed  Ion  milling 
where  the  sample  was  milled  at  room  temperature  until  perforation  occurred. 

In  liyht  of  the  ease  of  Introducing  artifacts  into  MCT  because  of  its  ease 
of  mechanical  deformation  even  at  room  temperature  [6],  It  is  to  be  emphasized 
here  that  this  sample  preparation  technique  was  attempted  because  of  Its  previous 
success  on  CdTe,  not  because  It  has  been  shown  to  be  successful  on  MCT. 

Indeed,  our  previous  experiences  w1tlf*1on  milling  showed  that  this  thinning 
technique  does  Introduce  defects. 

The  Ion  milling  employed  In  final  thinning  step  Is  Indeed  seen  to  Intro¬ 
duce  large  densities  of  defects.  The  results  of  this  preparation  technique 
are  shown  In  Figure  1.  Conventional  2  -  beam  Imaging  at  the  low  magnifi¬ 
cation  of  the  Interfacial  region  reveals  a  defect  density  that  Is  very  much 
higher  (2  to  3  orders  of  magnitude)  in  the  CdTe  than  In  the  MCT.  Figures  la,  b, 
and  c  are  micrographs  Illustrating  this  effect.  These  are  taken  from  the  interface 
region,  away  from  the  Interface  on  the  MCT  side,  and  away  from  the  interface 
on  the  CdTe  side,  respectively.  The  defects  In  the  CdTe  appear  to  consist 
of  dense  tangles  of  dislocations  and  large  numbers  of  small  loops.  The  defect 
density  decreases  over  several  microns  from  the  heterojunction  into  the  MCT 
where  the  defects  are  predominantly  single-line  dislocations  with  occasional 
large  loops.  Since  this  density  of  defects  is  extremely  high,  it  Is  seen 
that  this  current  preparation  technique  is  rendered  prohibitive  In  studying 
the  original  defect  structure  associated  with  the  hetero junction  since  it 
causes  damage  that  completely  masks  this  original  defect  structure. 


However,  since  the  MCT  and  the  CdTe  have  been  ion  milled  under  identical 
conditions  in  this  technique,  a  direct  comparison  of  the  ion  bombared  induced 
defects  between  the  MCT  and  CdTe  can  be  undertaken.  In  addition,  the  MCT-CdTe 
heterojunction  can  be  examined.  It  is  therefore  possible  to  extract  impor¬ 
tant  information  of  the  defect  energetics  from  these  ion-damaged  cross-sectional 
sampl es. 

Additional  insight  into  the  nature  of  these  ion-induced  defects  is 
obtained  if  higher  magnifications  are  used  in  studying  the  defects.  Lattice 
images  (to  our  knowledge,  these  are  the  first  such  images  obtained  from  HgCdTe) 
of  the  heterojunction  region  reveal  a  much  higher  density  of  dissociated  60° 
dislocations  and  interstitial  loops  in  the  CdTe  than  in  the  MCT  consistent 
with  the  above  observations  on  their  density.  This  high  density  is  commonly 
observed  in  ion-milled  CdTe  [7].  Typical  results  are  shown  in  Figure  2. 

Figure  2a  shows  a  lattice  image  of  a  representative  area  in  the  CdTe  side  of 
the  heterojunction  where  many  extrinsic  stacking  faults  (ESF)  are  readily 
observed.  That  there  exist  a  much  lower  density  of  defects  in  the  MCT  side  of 
the  heterojunction  can  be  seen  in  Figure  2b,  where  a  representative  area 
contains  no  defects  at  all.  Figure  2c  shows  an  almost  perfect  60°  in  the  MCT, 
where  the  arrow  is  pointing  down  the  extra  half-plane.  To  the  best  of  our 
knowledge,  perfect  60°  dislocation  have  not  been  observed  in  CdTe.  Figure  2d 
shows  a  typical  ESF  in  the  MCT,  with  the  width  between  the  two  Shockley  partial 
dislocations  always  observed  to  be  less  irv  this  material  than  in  CdTe. 

That  there  is  a  much  lower  defect  density  in  the  HgCdTe  than  in  the  CdTe 
when  damages  under  the  same  milling  conditions  at  first  seems  surprising. 

Cole  has  shown  that  HgCdTe  is  easily  deformable  even  at  room  temperature,  which 
produced  large  densities  of  dislocations  [6].  Cole  also  showed  that  the  rate 
of  glide  of  dislocations  is  similar  in  the  two  materials  because  the  room 
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temperature  hardness  of  CdTe  and  MCT  (x  *  0.3)  are  similar.  In  light  of 
this  information,  intuition  leads  us  to  expect  that  just  as  many  or  more 
defects  will  be  present  in  the  HCT  side  of  this  Ion-milled  heterojunction. 

That  this  Is  not  the  case  suggests  that  the  defect  energetics  are  different  for 
the  two  materials,  and  that  some  recovery  mechanism  exists  In  the  HCT,  the 
rate  of  which  being  much  more  rapid  In  this  material  than  in  CdTe.  The 
differences  in  the  defect  densities  on  either  side  of  the  heterojunction  cannot 
be  attributed  to  differences  in  the  glide  behavior  since  the  hardness  numbers 
are  similar,  therefore  other  recovery  mechanisms  must  be  occurring.  This  in 
agreement  with  theoretical  calculations  of  Harrison  [8]. 

Insight  Into  the  difference  of  recovery  rate  for  these  two  materials  can 
be  gained  from  these  results.  The  Information  obtained  from  the  lattice 
Images  taken  from  both  of  these  regions  reveal  that  there  Is  a  small  dissocia¬ 
tion  of  the  Shockley  partial  dislocations  In  the  HCT  relative  to  the  CdTe, 
which  suggests  that  the  climb  of  dislocations  In  HCT  Is  much  easier  than  In  the 
CdTe.  Dislocation  climb  Is  a  diffusion-controlled  process  with  vacancies  diffusing 
to  the  dislocation  core  as  the  rate  determining  step.  The  narrower  stacking 
fault  has  fewer  atomic  sites  to  which  vacancies  must  dlffust  to  cause  It  to 
climb,  hence  this  process  should  be  more  rapid  in  the  HCT  with  its  correspon¬ 
dingly  narrow  ESF.  This  might  also  be  aided  by  a  large  density  of  point 
defects  In  Ion  bombarded  HCT,  arising  from  the  weak  Hg-Te  bond  and  the  con¬ 
sequent  ease  of  formation  of  the  Hg  vacancy  [8,9]. 

Since  the  nature  of  the  damage  Introduced  by  ion  milling  is  similar  to 
that  of  Ion  Implantation  damage,  TEN  lattice  imaging  could  be  used  to  study 
the  nature  of  these  Ion-Implantation  defects.  Thus,  lattice  Imaging  might 
find  suitable  applications  In  understanding  some  of  the  critical  technolog¬ 
ical  problems  associated  with  Ion  Implantation  damage  to  the  HCT  lattice 
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as  well  as  those  Involved  in  LPE  growth  of  MCT  on  CdTe.  In  the  case  of 
studying  heterojunctions  in  cross  section,  methods  to  thin  the  material 
without  introducing  artifacts  must  be  developed  in  order  to  study  the  original 
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Conventional  2-beam  TEH  micrograph  of  the  ion-milled  HgCdTe-CdTe 
heterojunction  showing  the  defect  structure  of  this  damaged 
regi on . 

In  the  HgCdTe  away  from  the  heterojunction  showing  lower  density 
of  defects. 

In  the  CdTe  side  of  the  heterojunction  showing  2-3  orders  of 
magnitude  higher  density  of  defects. 

TEM  lattice  image  on  the  CdTe  side  of  the  heterojunction. 

This  representative  area  is  seen  to  possess  a  high  density  of 
extrinsic  stacking  faults  (ESF). 

Lattice  image  taken  from  a  representative  area  on  the  HgCdTe 
side  of  the  heterojunction  showing  no  defects. 

Lattic  image  showing  an  almost  perfect  60°  edge  dislocation  in 
the  HgCdTe  side  of  the  heterojunction.  The  arrow  points  down 
the  extra  half-plane. 

Dissociated  dislocation  in  the  HgCdTe.  Notice  the  much  narrower 
extrinsic  stacking  fault  (ESF)  width  in  this  material  as  com¬ 
pared  to  the  CdTe  ESF  width. 


Section  III-4 

Activated  Oxygen  Uptake  on  HgTe,  CdTe,  and  ^90.69^0.31 Ted 

J.  A.  Si1bermanb,  0.  Laserc,  I.  Lindau,  and  W.  E.  Spicerd 

Stanford  Electronics  Laboratory 
Stanford  University 
Stanford,  CA  94305 

J.  A.  Wilson 

Santa  Barbara  Research  Center,  Goleta,  CA  93017 

A  study  of  the  oxidation  of  the  semi  conductor  alloy  Hg]_xCdxTe  can 
afford  insights  into  the  surface  chemistry  of  the  alloy  and  the  way  in  which 
the  properties  of  the  constituents  (eg.  HgTe  and  CdTe)  are  modified  by 
alloy  formation.  From  experimental  Cl]  as  well  as  theoretical  [2]  inves¬ 
tigation  of  the  alloy  electronic  structure  it  is  known  that  the  cation 
s-electron  derived  contribution  to  the  deep  valence  states  in  Hgi_xCdxTe 
reflect  distinct  Hg  or  Cd  origin  rather  than  indicate  a  single  average 
or  virtual  cation.  The  separate  contribution  to  the  valence  electronic 
structure  by  the  cations  implies  a  difference  in  the  bonding  of  Hg  and 
Cd  to  Te.  The  more  metallic  nature  of  the  HgTe  bond  reduces  the  bond 
strength  in  HgTe  compared  to  CdTe  [3],  and  the  observation  of  binding 
energy  shifts  in  the  Cd  4d  and  Hg  5d  levels  in  Hgi_xCdxTe  compared  to 
the  binary  compounds  suggests  the  disparity  in  bond  strength  increases 
upon  alloy  formation  [1].  In  a  comparative  study  of  oxygen  uptake  on 
CdTe,  HgTe,  and  Hgo.69Cdo.31 Te  discussed  below,  the  modification  of  the 
Cd-Te  bonding  properties  are  manifest  in  a  decrease  in  activity  of  this 
site  in  the  alloy  compared  to  CdTe,  and  the  weaker  Hg-Te  bond  is  observed 
to  play  a  rate  determining  role  in  the  alloy  oxide  formation. 

To  examine  the  oxidation  chemistry  of  these  compounds,  photoemission 
spectroscopy  using  synchrotron  radiation  was  used  to  monitor  the  effects 
of  oxygen  exposure  on  the  surfaces  of  (110)  oriented  HgTe,  CdTe  and 


^90.69^0. 31  Te  cleaved  in  vacuum.  Because  these  compounds  are  highly 
inert  to  oxygen  in  the  molecular  ground  state  [4],  uptake  was  stimulated 
by  a  hot  filament  positioned  in  line-of-sight  of  the  sample  face.  So 
that  a  linear  exposure  scale  could  be  established,  each  sample  was  exposed 
to  multiples  of  a  fixed  dose  (5  x  10"®  torr  of  O2  for  500  sec)  thereby 
allowing  comparison  of  the  rate  of  oxidation.  The  photoemission  spectra 
of  the  valence  band,  Hg  5d,  and  Cd  4d  core  levels  (Figure  1)  and  the  Te  4d 
core  lines  (Figure  2)  following  each  exposure  document  the  interaction  of 
oxygen  activated  by  the  line-of-sight  hot  filament  with  Hg-|_xCdxTe.  Oxygen 
accumulates  with  each  exposure,  as  evidenced  by  the  increase  of  0  2p 
derived  emission  at  about  46  eV  (Figure  1).  A  second  oxygen  induced 
structure  occurs  near  42.2  eV  and  overlaps  the  Hg  5d5/2  line;  this 
second  feature  is  readily  observed  in  the  case  of  CdTe.  The  binding 
energy  of  the  main  0  2p  structure  relative  to  the  valence  band  maximum 
increases  by  roughly  0.5  eV  in  HgTe  compared  to  CdTe.  In  addition  to 
growth  of  the  oxygen  related  structure,  exposure  to  activated  oxygen 
results  in  attenuation  of  the  Hg  5d  core  levels.  (Because  of  the  substantial 
overlap  of  the  Hg  5ds/2  peak  with  the  oxygen  derived  structure,  this 
is  best  followed  with  the  Hg  5d3/2  line.)  The  Cd  emission,  on  the 
other  hand,  slightly  strengthens,  indicating  that  while  Hg  is  released 
upon  oxygen  uptake,  Cd  is  retained  in  the  oxide.  Release  of  Hg  is  observed 
in  the  initial  stages  of  HgTe  oxidation  as  well. 

Oxide  formation  is  signalled  from  the  initial  exposure  by  the 
appearance  of  Te  4d  emission  approximately  3.6  eV  more  tightly  bound 
-(Figure  2)  which  increases  in  magnitude  with  additional  exposure.  With 
each  dose,  the  Te  signal  at  the  clean  surface  position  is  attenuated  by 
the  growing  oxide  film.  That  the  Te  in  the  alloy  and  Hg  signals  diminish 


at  the  same  rate  suggests  Hg  Is  released  as  a  consequence  of  or  as  a  step  in 
the  oxidation  process  at  the  initial  stages  studied  here.  In  addition  to  the 
oxidized  and  lattice  bound  Te,  a  third  state  of  Te  results  from  the  first 
exposure.  Emission  from  4d  states  0.6  eV  more  tightly  bound  than  in  the 
bulk  produces  the  asymmetric  broadening  of  the  core  level  emission  near  the 
bulk  Te  position  seen  in  the  spectra  (Figure  2).  The  binding  energy  of  this 
third  contribution  to  the  Te  emission  indicates  elemental  Te  [6]  is 
produced  by  the  exposure  and  probably  results  from  a  loss  of  Hg  in  excess 
of  the  quantity  of  oxide  formed  during  the  first  dose.  The  extent  to 
which  this  occurs  varies  from  surface  to  surface  and  may  depend  on  the 
number  of  imperfections  produced  upon  cleavage. 

Examination  of  the  rate  of  oxygen  uptake  by  the  three  materials  may 
provide  further  information  as  to  the  reactions  controlling  oxide  growth. 

From  the  ratio  of  oxidized  to  unoxidized  Te  peaks  in  the  spectra  for  each 
sample  and  after  suitable  correction  for  the  inelastic  mean  free  path  of 
the  photoemitted  electron  (assumed  constant),  the  thickness  of  the  oxide 
film  can  be  determined  as  a  function  of  the  number  of  standard  exposures 
for  the  three  compounds.  This  data  is  given  in  Figure  3.  Several  obser¬ 
vations  regarding  the  relative  rates  relate  to  the  oxidation  process.  In 
every  case  the  oxide  grows  in  proportion  to  the  exposure.  In  the  avail¬ 
able  data,  there  is  no  evidence  of  saturation  for  films  grown  in  excess 
of  at  least  two  monolayers.  The  rate  of  oxide  production  for  the  alloy 
is  69%  as  great  as  that  of  the  binary  HgTe;  i.e.,  the  rate  of  growth 
scales  with  the  number  of  HgTe  bonds.  This  observation  is  consistent 
with  the  correlation  of  the  release  of  Hg  and  oxide  formation  noted  in 
the  core  level  spectra.  Finally,  that  the  alloy  of  this  composition 
oxidized  more  slowly  than  CdTe  implies  reduced  activity  at  the  Cd-Te 


sites  in  the  alloy.  The  rate  limiting  role  of  the  weaker  Hg-Te  bond  as 
well  as  the  modified  interactions  of  the  Cd-Te  sites  are  consistent  with 
the  changes  in  electronic  structure  between  HgTe,  CdTe,  and  the  ternary; 
studies  of  additional  compositions  as  well  as  identification  of  the 
nature  of  the  desorbing  Hg  will  provide  further  insights  into  the  oxidation 
reaction  and  thus  detail  the  relation  of  the  surface  electronic  structure 
and  oxygen  chemistry. 
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Figure  Captions 

Figure  1.  Photoemission  spectra  from  the  valence  band,  Hg  5d  (40.5  and 
42.3  eV),  and  Cd  4d  (centered  around  39.5  eV)  core  levels  for 
the  clean  surface  and  following  each  oxygen  dose.  The  valence 
band  maximum  is  located  at  50.2  eV. 

Figure  2.  Photoemission  spectra  of  the  Te  4d  core  lines  for  the  clean  surface 
and  with  increasing  oxygen  exposure. 

Figure  3.  Thickness  of  the  oxide  film  measured  in  multiples  of  the  electron 
Inelastic  mean  free  path  versus  exposure  in  multiples  of  a 
standard  dose. 
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Section  III-5 


Low  Temperature  Cleavage  of  MCT. 
A)  Motivation 


A  largo  body  of  work  ia  this  prograa  Am  boos  performed  on  tbo  (110) 
fact  of  llCT  foraod  by  cleaving  uaanaealed  oolid  state  recrystalizod 
material  ia  vacuum.  Doopito  tbo  1018  to  IO17  ca"3  p-type  character  of  the 
samples,  measurement  of  tbo  Fermi  lorol  pooitioa  at  tbo  oarfaco  Boot  oftoa 
iadicatoo  aa  a-typo  oarfaco.1  Tbo  typo  conversion  io  dao  to  defocto 
created  by  tbo  cleavage.  A  high  doaoity  of  dioloeatioao  has  beoa  observed 
at  tbo  odgoo  of  etepe  foraod  oa  tbo  cleavage  oarfaco  where  tbo  cleave 
jaaps  froa  oao  (110)  plaao  to  another.3  While  sitoo  asoociatod  with  sacb 
dioloeatioao  Bay  aot  tboaoolvoo  bo  electrically  active,  tbs  obsorvatioa  of 
theee  ortoadod  defocto  clearly  iadicatoo  tbo  oaorgy  to  oigaificaatly 
perturb  the  lattice  io  available  apoa  cleaving. 

Oao  aoaao  to  redace  tbo  possibility  of  defect  croatioa  is  to  redaco  the 
temperature  of  the  sample  at  the  time  of  cleavage.  Ia  a  study  of  the 
mechanical  properties  of  MCT3.  Cole  observed  a  strong  temperature 
dependence  of  the  yield  stress  for  tempo rataree  above  room  temperature. 
By  cooling,  the  material  may  harden  sufficiently  that  fracture  occurs  only 
along  a  single  (110)  plane. 

B)  Reealte  of  Initial  Attempt  at  Low  Temperature  Cleavage. 

> 

Aa  altrabigb  vacuum  ample  manipulator  for  use  with  the  angle-resolving 
chamber  was  fabricated  with  the  capability  to  cool  and  moderately  beat 
cleavable  samples  traasfered  into  the  vacuum  system  from  air.  So  that  tbs 
low  temperature  of  the  sample  could  be  maintained  at  the  time  of  cleavage, 
the  cleaver,  consisting  of  a  flat  annealed  copper  anvil  against  which  the 
sample  rests  and  a  sharp  tungsten  carbide  blade  driven  by  a  linear  motion 
feedthrough,  was  modified  so  that  liquid  nitrogen  cooling  could  be  applied 
to  the  anvil.  The  x»0.2  MCT  sample  used  achieved  a  minimun  temperature  of 
ISO  X.  Because  the  anvil  was  warmer,  the  temperature  at  the  time  of 
cleavage  was  *210  X.  The  resulting  surface  had  a  visible  morphology  and 
surface  potential  (the  surface  was  n-type)  which  fell  within  the 


i ,*'■ 

Era* 


distribution  of  results  achieved  to  dots  at  room  temperature. 

Tbs  constructed  manipulator  will  prows  of  grut  utility  in  extending 
studios  of  oxygon  uptaks  or  natal  dspooition  to  teaperaturee  otbsr  than 
room  tsnparaturs.  It  will  naks  possibls  sxpsriasnts  rslating  to  Rg  loss  or 
surfaes  structura  which  in  fact  raquira  hsating  or  cooling.  An  effort  will 
be  nads  to  iaprowa  the  psrfornancs  of  tha  cooling  systsa  of  the  cleaver 
and  to  repeat  the  low  temperature  cleaving  experiment. 
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Seotion  IV 
Alloy  Bonding 

Section  IV- 1 

A  number  of  calculations  of  properties  of  semiconductor  alloys  with 
emphasis  placed  on  the  II-VI  compounds  involving  HgTe  have  been  started. 

These  calculations  include: 

1.  determination  of  the  many  body  generalizations  of  the  atomic  term  values, 
electron  extraction  energies  and  many  body  excitation  energies  for  all  the 
group  II,  III,  IV,  V,  and  VI  elements  that  are  important  in  tetrahedral 
semiconductors.  Their  use  in  place  of  traditional  term  values  takes 
account  of  the  influence  of  some  important  correlation  effects  on  semi¬ 
conductor  bond  energies; 

2.  calculation  of  the  bond  energies  of  the  "pure"  semiconductor  compounds  and 
the  modifications  of  these  energies  in  pseudo-binary  alloys; 

3.  calculation  of  the  bond  length  changes  in  alloys  and  deviations  from 
Vegard's  law; 

4.  initiation  of  theory  of  vacancy  formation  that  includes  effects  due  to 
relaxation  and  charge  redistribution  around  the  impurity; 

5.  extension  of  Harrison's  theory  of  the  shear  coefficients  to  dislocation 
energies  in  pure  compounds  and  alloys; 

6.  devising  of  a  theory  of  the  hardness  of  semiconductors  and  their  alloys; 
and 

7.  working  out  an  improvement  in  the  theory  of  the  energy  positions  (deep 
states)  of  impurities  in  pure  semiconductors  and  alloys. 

Not  all  these  calculations  are  completed,  but  several  are  in  a  stage  worth 
reporting.  In  Section  IV-2  the  current  status  of  our  work  on  bond  energies 
and  lattice  constants  will  be  presented.  Section  IV-3  has  a  discussion  of 
dislocation  energies  and  hardness.  Conclusions,  recommendations  based  on 
them,  and  future  plans  are  presented  in  Section  IV-4. 
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Section  IV-2 


"*  BOND  ENERGIES  AND  LENGTHS  IN  AN  ALLOY 

The  theory  of  the  bond  energy  shifts  in  an  alloy  that  was  presented  in 
our  mid-year  report ^  has  now  been  refined  and  extended.  All  the  important 
III  -  V  and  II -VI  compounds  have  been  treated.  The  details  of  these 
calculations  will  be  prepared  as  a  technical  paper.  The  most  important 
results  are  summerized  in  Figures  (II  -  1)  and  (II  -  2)  where  the  shift  of  a 
given  type  bond  energy  (e.g.  HgTe)  due  to  having  one  bond  of  a  different 
type  (e.g.  CdTe)  in  one  of  the  three  adjacent  bonds  are  presented.  These 
energy  shifts  for  now  Include  only  the  largest  terms;  those  that  arise 
from  local  charge  redistribution  in  the  alloy.  There  are  smaller  terms  due 
to  bond  length  and  bond  angle  distortions  that  must  still  be  included  to 
complete  this  calculation. 

Figure  II  -  1  is  a  compilation  of  these  bond  energy  shifts  for  the  III-V 
compounds  and  Figure  II  -  2  is  a  similar  plot  for  the  II  -  VI  compounds.  For 
each  bond  type*  e.g.  GaAs,  the  lines  under  a  cation  denote  the  energy  shift 
if  one  of  the  three  neighboring  bonds  coupled  to  the  anion  have  its  cation 
replaced  by  different  atom  e.g.  Al.  A  positive  energy  shift  is  destablizing 
and  a  negative  shift  is  stablizlng.  An  ordering  Al*Ga*In  for  the  cations* 
and  P,As»Sb  for  the  anions  is  used.  An  arrow  under  Ga  pointing  to  the 
right  means  In  and  to  the  left  means  Al.  A  double  arrow  means  move  two 
steps  In  the  Indicated  direction.  Thus  in  the  alloy  Al^_xGaxAs  according 
to  this  figure  the  GaAs  bond  is  destablized  by  ~43meV  by  an  AlAs  neighboring 
bond*  and  the  AlAs  bond  is  stablized  by  ~46meV.  If  there  are  two  neighboring 


bonds  with  the  different  atons  these  numbers  roughly  double,  and  for  three 
they  triple. 

Focus  now  on  Fig.  (II-2_  and  In  particular  the  HgTe  bond.  This  bond 
Is  destabllzed  by  ~  30  meV  per  neighboring  CdTe  bond  and  ~5  meV  per  neighboring 
HgS  or  HgSe  bond.  However,  the  bond  Is  stabllzed  by  ~  16  meV  per  neighboring 
ZnTe  bond.  Since  according  to  these  calculations  the  bond  strength  of  HgTe 
is  only  -0.480  eV,  In  a  Hg^Cd^Te  alloy  where  a  given  HgTe  bond  is  surrounded 
by  three  CdTe  bonds  its  bond  energy  is  reduced  by  25Z.  This  can  have  a 
substantial  effect  on  a  Hg  vacancy  formation  energy.  We  have  estimated  the 
Hg  Schottky  defect  formation  energy  In  HgTe  to  be  0.78  eV,  and  the  presence 
of  Cd  could  reduce  this  to  as  little  as  0.3 eV.  It  is  unlikely  to  find  this 
large  a  reduction  in  an  alloy  with  x  a  0.2.  However,  in  an  alloy  with  x  near 
unity  the  Hg  vacancy  formation  energy  could  be  quite  small.  This  could  have 
dlsasterous  consequences  to  the  stability  of  CdTe  solar  cells  where  a  low 
concentration  of  HgTe  is  added  to  slightly  reduce  the  bandgap.  While  the 
consequences  are  not  so  severe  In  the  concentration  range  used  in  Infrared 
detector  applications,  they  are  still  potentially  troublesome.  We  suspect 
this  effect  helps  account  for  the  ease  with  which  Hg  vacancies  are  found 
to  form  in  the  alloy. 

If  one  calculates  the  bond  length  changes  due  only  to  effects  caused 
by  the  charge  shifts  which  are  driven  principally  by  the  differences  in  the 
ionic  properties  of  the  bonds  then  they  all  shift  in  the  wrong  direction 
compared  to  experiment.  Experiments  on  GalnAs  have  shown  that  in  the  alloy 
the  bond  lengths  of  the  pure  materials  tend  to  be  preserved  but  the  shorter 


GaAs  bond  expands  slightly  (~  21)  and  the  longer  InAs  bond  contracts  by 
~2Z.  The  charge  shift  contribution  yields  the  correct  order  of  magnitude 
but  the  wrong  sign.  The  charge  shift  term  has  a  large  effect  on  the 
metallic  contribution  to  the  bond  energy.  But  this  contribution  is  a 
small  one.  There  is  Inherent  in  the  bond  length  change  the  necessity  for 
bond  angle  distortions  and  a  change  in  the  average  bond  lengths  to  allow 
the  lattice  to  fit  together.  These  bond  angle  distortions  cause  small 
changes  in  the  covalent  contribution  to  the  bond  energy,  but  this  contri¬ 
bution  is  a  large  one.  Once  the  bond  angle  distortion  terms  are  included 
the  bond  length  changes  tend  to  reverse  their  sign.  To  obtain  answers 
the  calculation  must  be  done  self  consistently,  permitting  the  average 
lattice  spacing  to  take  on  its  proper  value.  Instead  of  simple  algebraic 
expressions,  the  answers  now  require  a  computer  program  that  itterates 
to  self  consistent  results.  We  shall  report  the  details  of  these  calculations 
when  they  are  completed. 

If  one  adopts  a  model  in  which  the  spring  constants  associated  with 
a  bond  are  supposed  to  retain  their  pure  crystal  values  in  an  alloy,  then 
it  can  be  demonstrated  that  the  bond  length  shifts  have  the  observed  sign 
and  magnitude.  This  model  also  accounts  for  the  observed  small  deviations 
from  Vegard's  law  in  which  the  average  bond  length  as  a  function  of  alloy 
concentration  bows  down  slightly  from  a  straight  line  projection.  At 
present  we  have  no  theoretical  justification  for  a  conserved  spring  constant 
model,  but  the  fact  that  it  gives  the  correct  answers  is  instructive. 


Section  IV- 3 


SEMICONDUCTORS  AND  THEIR  ALLOYS: 
DISLOCATIONS  AND  HARDNESS 


Dislocations  in  semiconductors  are  detrimental  to  device  function, 
serve  as  channels  for  impurity  migration  and  trapping  which  causes  non  uni¬ 
form  doping  and  degrades  p-n  junctions.  They  also  decrease  the  material's 
resistance  to  plastic  deformation.  The  aim  of  this  paper  is  to  provide 
insights  into  the  underlying  physical  mechanisms  controlling  dislocations 
and  the  hardness  of  semiconductors  and  their  alloys;  and  then  to  suggest 
strategies  for  decreasing  dislocation  densities.  It  is  well  established 


that  the  hardness  of  tetrahedrally  coordinated  semiconductor  materials. 

Group  IV,  III-V  and  II-VI  compounds  exhibit  a  sharp  variation  with  their 

_9 

near  neighbor  distance  d,  approximately  proportional  to  d  for  one  group 
(2) 

of  seven  compounds.  7  (The  lattice  constant,  a,  and  d  are  related  by 
d  ■  /Ja/4.)  Thus,  semiconductors  with  small  lattice  constants  tend  to  be 

(3) 

harder  materials.  These  same  materials  have  larger  stiffness  coefficients 
and  have  fewer  dislocations  in  as-grown  crystals.  The  composition 


x  variation  of  the  hardness  of  semiconductor  alloys  A^_xBx  bows  up  from  a 

straight  line  projection  between  the  respeceive  hardnesses  of  the  pure 
(2) 

constituents;  The  composition  variation  of  the  lattice  constant  of  alloys 


bows  down  from  the  Vegard's  rule  straight  line  projection  between  those  of 
the  constituents.^  we  shall  now  interrelate  these  observations  and  some 


other  more  recent  measurements,  through  an  extension  of. Harrison's  theory  of 

(8) 

the  bonding  of  these  semiconductors.  ’ 


The  major  reason  for  the  properties  just  described  can  be  understood 
qualitatively  as  follows.  These  materials  have  a  low  coordination  number 


*  P 


with  only  four  near  neighbors.  There  are  three  contributions  to  their  bonds. 

The  most  Important  one  is  the  covalent  energy  that  results  from  the  inter- 

3 

actions  between  the  sp  hybrids  on  the  different  sites.  The  strength  of 
this  Interaction  is  directional  dependent  and  bond  angle  distortions  require 
energy.  The  ionic  contribution  to  the  bond  strength  is  fairly  isotropic  and 
adds  little  to  the  shear  coefficients  per  bond.  In  crystals  with  higher 
coordination  the  ionic  contribution  to  shear  coefficients  is  more  important. 
The  metallic  contribution  to  the  bond  energy  is  in  this  case  also  bond  angle 
dependent.  In  metals  (usually  fee,  or  bcc  structures)  with  high  coordination 
numbers  this  term  dominates  the  shear  coefficients.  The  shear  coefficients 

(combinations  of  and  C^-C^)  depend  on  crystal  orientation  and  in 

(T\  3*2  2  -2 

Harrison's  notation'  1  are  proportional  to  +  V3  )  and  V2~  d  . 

V2  is  the  covalent  and  V3  the  ionic  energy  in  this  notation.  Thus  in  a  case 
where  the  covalent  energy  dominates  (V2  »  V^)  the  shear  coefficients  are 
proportional  to  d”^,  while  in  the  other  limit  they  are  smaller  and  propor¬ 
tional  to  d”^  . 


There  are  two  contributions  to  the  dislocation  energy  per  unit  length v 

a  short  ranged  core  energy  and  a  long  range  elastic  energy  arising  from  an 

Integral  over  the  energy  of  small  bond  angle  distortions  distributed  over 
(9) 

many  sites.  The  elastic  contribution  is  typically  ~  5  times  larger  than 
the  core  energy  and  while  the  elastic  term  dominates,  the  core  should  not  be 
completely  neglected.  However,  in  this  treatment  where  the  major  trends  are 
being  sought  we  shall  do  so.  The  distortion  energy  per  unit  length  is  pro¬ 
portional  to  the  product  of  a  shear  coefficient  and  the  square  of  the  Burger's 
(9) 

vector.  The  magnitude  of  the  Burger's  vector  is  porportional  to  the  lattice 

3 /2 

constant,  hence  ~  +  V^) 
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Hardness  Is  determined  by  applying  a  known  force,  F,  to  a  probe  of  a 
prescribed  shape  driving  it  into  the  surface  of  the  sample. ^  The  area.  A, 
of  the  resulting  Indentation  is  measured,  and  the  hardness  is  the  force  per 
unit  Indented  area.  Many  dislocations  must  be  formed  to  allow  the  probe  to 
indent  the  semiconductor.  If  the  indenter  is  a  rectangular  pyramid  where 
the  angle  between  the  normal  to  the  point  and  an  edge  is  6,  then  the  major 
contribution  to  the  hardness  is 


F  _  Fh  _  fT 
A  Ah  Ah 


(3-1) 


where  c.j.  is  the  work  required  to  cause  the  Indenter  to  penetrate  to  a  depth  h. 

A  side  view  of  the  Indentation  in  a  cut  through  its  center  is  Illustrated 
schematically  in  Figure  (3-1).  The  top  of  the  indentation  has  side  length  W 
so  A  -  W^.  The  Burger's  vector  has  magnitude  b  which  is  approximately  equal 
to  the  bond  length  d.  The  number  of  dislocations  required  to  accommodate 

an  indentation  to  depth  h  is 

«  _  h  _  W  cote  (3-2) 

Nh  b  ~"""2b 

where  e  is  the  angle  between  the  normal  to  the  tip  of  the  Indenter  and  a  side. 
The  edges  of  the  extra  atom  planes  that  are  driven  from  the  Indented  volume 
into  the  bulk  of  the  semiconductor  are  shown  as  lines  in  the  figure.  The 
ends  of  these  lines  are  the  positions  of  the  dislocations  which  are  per¬ 
pendicular  to  the  plane  of  the  figure.  The  planes  driven  to  the  sides  each 
have  a  finite  extent  and  a  trapazoidal  shape.  The  planes  dlrven  down  under 
the  Indenter  have  a  square  shape. 

Most  of  the  work  done  on  the  indentor  goes  into  the  energy  to  form  the 
dislocations.  There  are  two  major  contributions  to  this  energy.  The  first 
is  the  energy  needed  to  form  each  dislocation  and  the  second  is  the  interaction 
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energy  among  these  dislocations  .  Since  the  extra  planes  driven  to  the 
sides  of  the  indentation  have  a  finite  extent  there  are  both  edge  dislocations 
at  their  base  and  screw  dislocations  associated  with  their  termination.  The 
square  planes  driven  below  the  indentation  have  edge  dislocations  around  the 
sides  and  screw  dislocations  at  the  corners  to  make  the  turns.  Moreover, 
there  are  Interactions  among  the  dislocations  which  can  produce  a  minimum 
energy  configuration.  For  the  kind  of  situation  depicted  in  Figure  (3-1) 
it  always  costs  energy  to  position  a  second  dislocation  on  a  parallel  glide  plane 
to  one  already  present.  However,  the  magnitude  of  this  extra  energy  can  be 
minimized,  and  for  proper  configurations  of  the  dislocations  there  are 
attractive  forces  along  the  glide  planes  that  will  tend  to  position  the 
dislocations  into  the  minimum  energy  configuration.  The  minimum  configuration 
arises  when  half  the  atoms  from  the  Indented  volume  go  respectively  to  the 
side  and  below  the  lndenter.  Then  in  both  regions  the  maximum  angle  made 
between  successive  close-spaced  dislocation  lines  and  their  glide  planes  is 
cp  ■  tt/2  -  e/2  as  shown  in  the  figure.  This  is  the  minimum  energy  configuration. 
As  we  shall  demonstrate  shortly  the  minimum  energy  configuration  is  evidently 
established  when  soft  materials  are  Indented,  but  hard  ones  are  more  apt  to 
align  in  a  configuration  closer  to  a  maximum  energy  configuration  where  all 
the  atoms  go  either  to  the  sides  or  below  the  lndenter. 


An  approximate  expression  for  the  energy  required  to  indent  the  material 
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In  the  first  term,  the  energy  per  unit  length  to  form  an  eige  dislocation  is 


Gb 


Ei  "  A'^l-vV  WV 


(3-4) 


the  shear  coefficient  is  G,  the  Burger's  vector  b  *  d//3  for  an  indentation 
along  a  (100)  axis,  the  range  of  the  dislocation  R  is  taken  equal  to  W  (for 
want  of  a  better  approximation)  rQ  d  is  the  dislocation  core  radius,  the 
Poisson  ratio  v  ~  0.2  for  most  semiconductors,  and  is  the  length  of  the  1' 
dislocation.  In  the  second  term,  is  the  Interaction  energy  per  unit 
length  between  dislocations  i  and  j  assuming  they  have  parallel  glide  planes 
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r? is  the  separation  between  dislocation  1  and  j ,  and  cp^j  is  the  angle  a 
line  perpendicular  to  the  dislocations  and  joining  them  makes  with  the  glide 
plane  as  shown  in  Fig. (3-1).  Since  the  various  dislocations  in  a  region 
have  different  lengths,  the  net  interaction  energy  is  approximated  by  multi¬ 
plying  the  energy  per  unit  length  by  the  length  of  the  shorter  one.  The  1/2 
multiplying  this  term  prevents  double  counting  and  the  prime  on  the  sum 
indicates  that  terms  with  i  ■  j  are  to  be  omitted.  The  upper  limits  on  the 
sums  N  are  the  number  of  dislocations  in  one  region  (side  or  bottom)  associated 
with  one  edge.  For  the  minimum  energy  configuration  N  -  N^/2 .  The  four  that 
multiplies  the  bracket  accounts  for  the  four  sides,  and  the  two  for  the  two 
regions  for  each  side. 


We  now  encounter  our  first  surprise.  The  energies  per  unit  lentgh 

2 

and  have  comparable  magnitudes.  Since  there  are  ~  N  terms  in  the  inter¬ 
action  energy  sum,  only  N  terms  in  the  formation  energy  sum,  and  N  is  a  large 
number  typically  ~  10^,  the  interaction  energy  completely  dominates  the  hardness. 
Hence  neglecting  terms  due  to  screw  dislocations  core  energies,  and  other 


effects  that  lnfluanca  the  first  tars  are  unimportant.  However,  care  must 
be  taken  with  the  Interaction  terms.  In  Eq.  3-3,  a  number  of  secondary 
interactions  are  neglected,  some  positive  and  others  negative.  These  will 
he  added  later.  The  principal  neglected  terms  are  the  Interactions  between  two 
regions  on  each  side  (positive),  and  the  Interactions  between  adjacent  side  (and 
bottom)  regions  (positive,  and  the  Interaction  between  the  opposite  sides  In  the 
bottom  regions  (negative).  Comparison  of  the  results  with  experiment  will 
tell  us  how  Important  these  neglected  terms  are  likely  to  be. 


The  length  of  the  i'th  dislocation  will  be  taken  as 


Li  “  w  h  *  *  "  • 


(3-6) 


This  is  the  largest  length  of  the  side  inserted  planes  and  its  choice  partially 
accounts  for  interactions  between  the  otherwise  neglected  screw  dislocations. 


The  distance  rji  *  ri  ”  is  given  by 


r  -  V2b|i  -  Jj  ;  i,J  -  1,2,3 . N 


(3-7) 
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for  the  minimum  energy  configuration  and  a  tetrahedrally  bonded  semiconductor. 


Finally  in  the  indicated  configuration 


0 

■  9  ■  ^  -  2  for  all  i,  and  j. 

Inserting  these  expressions,  Eqs.  (3-5),  (3-6),  (3-7),  and  (3-8), 

into  the  sums  in  Eq.  (3-9)  yields  expressions  of  the  form 
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where 


g  e  V2b/R  «  1 


(3-9) 
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and  only  terms  of  order  N*  are  retained  in  Eq.  (3-9). 
Similarly  one  finds 

r  i-i  n 

2}ln  {lj.Lj}  -  |  22  OH)  +  (N-1)  22  1 

j  1  j-i+1 


S'  I  N2  -  Ni  +  l2/2 


(3-11 


Combining  these  results  yields 
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and  for  the  minimum  hardness 


a  .  G  cot 
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One  can  also  get  a  number  for  the  hardness  of  a  dislocation  con¬ 
figuration  in  which  all  the  material  is  either  pushed  to  the  side  or  the 
bottom.  In  this  case  the  factor  of  2  in  front  of  Eq.  (3-3)  is  removed, 

N  -  N.  ,  g  *  y2b/R,  and  cp  ■  tt/2  -  0.  Then  a  non-equilibrium  maximum  hardness 


in  the  context  of  this  model, is  obtained  (remember  that  some  terms  have  been 


neglected) 


G  cot  6  T  „  /cot  e^  ■ 

PH -75-/  +  12 


+  2  sin  e 


(3-14) 


The  proper  answer  for  most  materials  lies  somewhere  between  and 

For  an  lndenter  with  0  ■  tt/4  we  have 


H  .  -  0.0414  ~ 

min  1— v 


(3-15) 


and  H  /H  .  -  5.80.  Harrison  '  has  shown  that 

max  min 
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where  -fi  is  Planck's  constant  over  2tt,  m  is  the  free  electron  mass 

2  2  1/2 

(fi2/m  -  1.220  x  10-38  [joules  -m2]),  <*c  is  the  covalence  orc  -  V2/(V2  +  V3  ) 
and  d  is  the  bond  length.  Combining  Eqs.  (3-14)  and  (3-16),  using  v  ■  0.2, 
and  changing  the  dimensions  to  those  in  terms  of  which  experimental  hardness 
numbers  are  customarily  quoted  gives 


1.23  X  10 
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(d  [1]  )5  j^(mm)2  J 


(3-17) 


Values  calculated  from  Eq.  (3-17)  and  H  along  with  the  corresponding 

max 

experimental  values  Hgxp  are  given  in  Table  (3-1)  and  the  theoretical  versus 
the  experimental  results  are  plotted  in  Fig.  (3—2)  for  a  number  of  typical 
semiconductors . 


The  plot  in  Fig.  (3-2)  has  the  minimum  and  maximum  theoretical  values 
of  H  at  the  two  ends  of  the  vertical  lines  for  each  compound  plotted  as  a 
function  of  the  corresponding  experimental  values.  If  the  theory  were  perfect 
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Experimental  and  Theoretical  Values  of  the  Hardness  of  Several  Semiconductors  in  Units  of 
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and  the  experimental  values  were  accurate  the  points  would  £all  on  the 

indicated  unity  slope  line.  There  are  several  conclusions  that  can  be 

drawn  from  Fig.  (3-2).  First  the  order  of  magnitude  of  the  predicted  and 

measured  values  are  the  same.  Second  the  trends  from  one  compound  to  another 

are  properly  given  by  the  theory.  The  H  ^  values  are  generally  too  small, 

but  they  fit  the  soft  materials  better,  and  H  fits  the  harder  materials 

max 

better.  As  we  have  pointed  out  some  terms  have  been  neglected  that  will 

very  likely  increase  H  ,  (and  H  ),but  they  should  be  smaller  than  those 

mm  max 

retained. 


The  semiconductor  alloys  exhibit  a  composition  variation  of  their  hard- 

(2) 

ness  that  bows  up  relative  to  the  straight  line  virtual  crystal  behavior. 
Moreover  the  peak  in  the  bowed  curve  is  in  general  shifted  toward  the  higher 
hardness  side  of  the  curve.  The  bond  length  of  alloys  had,  until  recently, 
been  thought  to  vary  nearly  linearly  with  composition  (Vegard’s  law)  between 
those  of  the  alloy  constituents.  The  deviations  of  the  average  bond  length 
from  the  simple  concentration  weighted  average  behavior  has  a  downward 
curvature.  However,  new  EXAPS  measurements  on  Ga^  ^In^As .clearly  display  an 
entirely  new  behavior.  The  bond  lengths  of  the  constituents  remain  nearly 
constant  in  the  alloy.  While  the  longer  one,  GaAs,  becomes  shorter  and  the 
shorter  one,  InAs,  becomes  longer  as  the  alloy  composition  changes,  they  each 
change  by  less  than  2Z  of  their  pure  crystal  values.  In  order  for  the 
crystal  to  accommodate  the  resulting  microstrains  and  for  an  alloy  in  which 
the  cation  is  the  substitutional  constituent,  the  cation  sublattice  retains 

its  face  centered  crystal  structure  and  the  anions  adjust  their  positions 
according  to  their  near  neighbor  configuration.  For  example,  if  an  As  cation 

has  three  In  neighbors  and  one  Ga  neighbor,  it  will  shift  along  a  (1,1,1)  direc¬ 
tion  toward  the  Ga  neighbor.  One  can  shov  that  for  an  A^  alloy  the  average 

bond  length  3  is  given  to  lowest  order  by 

90 


d  -  a-*)  dA  +  X  dB  +  3x  (1-x)  (Ad*  +  AdJ)  (3-18) 

where  dj  Is  the  bond  length  of  the  pure  compound  of  type  j  ■  A,B,  and 
A4^  is  the  shift  in  bond  length  of  bond  A  in  a  configuration  in  which 
there  are  three  A  bonds  and  one  B  bond  (an  A(3)B(1)  configuration).  The 
first  two  terms  are  Vegsrd's  law  which  says  that  the  bond  lengths  vary 
linearly  between  those  of  the  pure  constituents.  The  correction  bowing 
term  depends  on  the  algebraic  sum  of  the  bond  length  changes  for  A(3)B(1) 
and  A(1)B(3)  configurations.  We  have  done  a  literature  search  on  measure¬ 
ments  of  d.(x)  and  found  that  in  all  cases  where  deviations  from  Vegard's  law 
have  been  determined  the  average  lattice  constant  bows  down,  so  a<^  +  Ad^ 
is  negative.  This  correlates  well  with  the  experimental  finding  that 
hardness  bows  upward  since  from  experiment  H  is  roughly  proportional  to 

-9 

d  .  Thus  a  large  part  of  the  bowing  of  the  hardness  of  semiconductor 
alloys  (sometimes  called  alloy  hardening)  correlates  with  the  extra  energy 
required  to  form  dislocations  and  other  mechanisms  e.g.  dislocation  pinning 
at  Impurity  sites,  often  Invoked  to  explain  this  kind  of  phenomena  In 
metals  probably  play  a  minor  roll  in  semiconductors. 

Dislocations  are  often  found  in  materials  as  they  are  grown.  Their 
density  is  determined  by  the  strains  to  which  they  are  subjected  In  the 
growth  process,  and  the  cool  down  schedule.  A  dislocation  constitutes  an 
excitation  relative  to  the  perfect  crystal  ground  state.  At  the  elevated 
grown  temperature  the  number  of  dislocation  present  is  controlled  by 
equilibrium  thermodynamics.  As  the  crystal  is  cooled,  a  critical  tempera¬ 
ture  is  reached  where  the  time  required  for  the  dislocations  to  equilibrate 
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becomes  long  compared  to  annealing  times  so  the  population  of  dislocations 
is  frozen.  For  some  systems  the  critical  temperature  may  be  the  growth 
temperature.  In  any  case  one  expects  materials  with  similar  growth  tempera¬ 
tures  to  have  fewer  dislocations  if  is  larger.  This  expectation  is 
confirmed  in  the  recent  work  reported  on  Zn^^Cd^e  bulk  material.  The 

best  CdTe  that  has  even  been  grown  has  dislocation  densities  in  excess  of 
5  -2 

5  x  10  cm  .  The  addition  of  only  4Z  Zn  reduced  the  dislocation  count  to 
less  than  5  x  104  cm”  .  The  ZnTe  bond  length  is  2.643  A  while  the  CdTe  bond 
length  is  2.805  A,  a  6Z  difference.  This  6Z  difference  in  bond  length 
translates  into  a  2Z  difference  in  the  dislocation  energy  per  unit  length 
for  1-x  «  0.04.  Dislocation  energies  for  typical  lengths  are  10  eV,  so 
a  2%  Increase  corresponds  to  ~  0.2  eV.  Using  a  critical  dislocation  freeze 
in  temperature  of  800°C  this  yields  a  population  decrease  of  ~  20.  These 
estimates  demonstrate  that  the  observed  factor  of  10  reduction  in  the  dis¬ 
location  density  of  the  Zr^  q£&q  g^Te  alloy  relative  to  that  of  CdTe  lies 
in  the  range  of  values  that  can  reasonably  be  expected. 

The  argument  just  presented  naturally  leads  to  a  strategy  for  decreasing 
dislocations  in  other  semiconductors.  If  an  alloy  is  made  of  the  material 
of  interest  with  another  compound  with  a  shorter  bond  length  then  the 
dislocation  density  should  be  reduced.  For  example,  this  suggests  that  the 
addition  of  a  small  amount  of  GaP  (d  -  2.359  A)  may  significantly  reduce  the 
dislocation  density  of  bulk  grown  GaAs  (d  ■  2.448  A).  It  has  been  demon¬ 
strated  that  the  addition  of  ~  1%  GaN  (d  ■  1.946  A)  to  GaAs  can  yield  a  larse 
volume  of  dislocation-free  material. 

We  have  demonstrated  that  the  dislocation  energies,  and  hardness  of 
tetrahedrally  bonded  semiconductors  are  rapid  functions  of  the  reciprocal 
of  the  bond  length.  The  bond  lengths  of  the  different  bonds  in  semicon- 


due Cor  alloys  nearly  retain  Chair  pure  eyrsCal  values,  and  the  average 
bond  length  always  bows  down  from  Vegard'a  rule.  This  observation  correlates 
well  with  the  finding  that  the  composition  variations  of  the  hardness  of 
semiconductor  alloys  bows  upward.  Moreover,  the  rapid  d  dependence  of 
dislocation  energies  explains  the  dramatic  decrease  of  the  dislocation 
density  in  bulk  grown  Zn^  g^CdQ  ^Te  material  relative  to  that  found  in  CdTe. 


Section  IV- 4 


IV  CONCLUSIONS  AND  FUTURE  PLANS 

A  theory  of  a  number  of  material  structural  properties  of  semi¬ 
conductor  alloys  has  been  established.  This  theory  accurately  predicts  many 
of  the  observed  trends,  and  even  some  quantitative  results.  For  the  focus 
material  of  this  study  Hg^^Cd^Te  our  major  conclusions  are 

•  Cd  destablizes  the  already  weak  HgTe  bond  and 
has  little  effect  on  the  dislocation  energy. 

e  Zn  both  stabilizes  the  HgTe  bond  and  increases 
the  dislocation  energy. 

More  generally  for  other  alloys  an  Important  result  is  that  if  one  has  a 
material,  e.g.,  GaAs,  and  one  wishes  to  reduce  the  dislocation  density  is 
to  alloy  it  with  another  material,  e.g.  GaP,  GaN,  BAs,  with  a  shorter  bond 
length.  There  is  experimental  evidence  for  several  alloy  combinations  that 
confirms  the  utillzy  of  this  strategy^ Obviously  one  must  not  go 
to  such  a  high  alloy  concentration  that  the  electronic  properties  of  the  host 
are  badly  degraded. 

The  plans  for  our  future  work  Include  the  following  casks: 

1)  Complete  the  alloy  bond  length  shift  calculation, 

2)  Complete  the  vacancy  formation  energy  calculation 
and  extend  it  to  find  diffusion  coefficients, 

3)  Improve  the  accuracy  of  the  hardness  theory  and 
extend  it  to  alloys, 

4)  Determine  the  mechanism  driving  the  thick  region 

of  misfit  dislocations  found  between  CdTe  substrates 
and  epitaxial  HgCdTe  layers. 


Establish  the  relation  between  alloy  bonding  properties 
and  optical  phonon  mode  frequencies. 

Begin  calculating  the  structural  properties  of  alloy 
surfaces. 

Finish  calculating  the  energy  levels  of  deep  impurity 
states  in  HgCdTe  alloys. 
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